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IS.  abstract 


The  present  theory  for  predicting  the  dependence  of  the  force 
of  adhesion  on  zeta  potential,  surface  charge,  electrolyte  or  sur¬ 
factant  concentration  and  pH  is  reviewed.  It  is  shown  to  be 
inadequate  for  explaining  experimental  observations.  New  theory  is 
developed  which  explains  these  observations  and  is  in  good  quan¬ 
titative  agreement  with  experimental  data  in  the  literature. 
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ABSTRACT 


The  present  theory  for  predicting  the  dependence  of  the  zeta 
potential  and  surface  charge  of  colloidal  particles  on  the  con-  . 
centration  of  electrolytes,  surfactants,  and  pH  is  reviewed  and 
shown  to  be  inadequate  for  explaining  experimental  observations. 

A  new  theory  is  developed  which  explains  these  observations  and 
is  in  gocd  quantitative  agreement  with  experimental  data  in  the 
literature . 

This  theory  considers  the  various  ways  in  which  a  surface 
charge  can  develop  on  a  particle.  These  include  adsorption  of 
complex  ions  derived  from  electrolytes  and  from  the  particle 
itself  and  dissociation  of  surface  groups  due  to  chemical  reac¬ 
tions  at  the  particle  surface.  Adsorption  isotherms  are  developed 
for  each  case.  Since  the  surface  charge  and  potential  enter  into 
the  expressions  for  these  isotherms,  the  relationship  between  the 
zeta  potential,  surface  potential  and  surface  charge  is  inves¬ 
tigated  .v  The  result  is  an  analytical  expression  for  planar  sur¬ 
faces  and  a  computer  program  for  spherical  surfaces."  When  the 
adsorption  isotherms  are  used  in  combination  with  these  results, 
the  manner  in  which  the  zeta  potential  and  surface  charge  depend 
on  ionic  strength,  surfactant  concentration  and  pH  is  completely 
described. 

The  theory  is  compared  with  experimental  data  available  in 
the  literature  and  good  agreement  between  them  is  found.  A  com¬ 
puter  program  was  developed  to  perform  the  regression  analysis 
required  for  estimating  the  adsorption  energies  and  number  of 
adsorption  sites  which  are  the  parameters  in  the  theory. 

The  theory  for  the  interaction  energy  between  two  particles 
is  also  investigated.  It  is  concluded,  that  the  theory  for  van 
der  Waala -London  forces  is  sufficiently  complete  so  that  no 
improvements  are  needed.  However,  in  case  of  the  electrical 
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forces  arising  from  double  layer  interaction,  the  present  theory 
is  extended  to  obtain  approximate  analytical  expressions  which 
give  this  interaction  energy  for  arbitrary  surface  potentials, 
particle  size,  and  separation  distance.  Approximations  available 
in  the  literature  are  restricted  to  either  particles  at  the  same 
potential  or  to  potentials  smaller  than  50  mv.  Although  exact 
numerical  calculations  have  been  performed,  the  results  are  in 
tabular  form  and  are  inconvenient  to  use.  Since  the  interaction 
energy  curve  is  to  be  used  Inter  in  the  program  in  further 
integrations,  it  is  necessary  to  have  analytical  expressions  to 
avoid  unwielding  numerical  problems. 
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SUBMICRON  SEPARATION  AND  DATA 


1.  INTRODUCTION 

As  motivation  for  why  the  work  to  be  presented  was  performed, 
two  unexplained  experimental  observations  which  are  commonly 
encountered  in  the  study  of  the  stability  of  colloidal  systems 
will  be  discussed.  Then,  the  present  theory  available  for  inter¬ 
preting  colloidal  phenomena  will  be  reviewed,  and  it  will  be 
shown  that  it  is  unable  to  account  for  these  experimental  results. 

The  most  typical  kind  of  data  that  can  be  obtained  for 
characterizing  a  colloidal  system  is  the  zeta  potential  versus 
concentration  curve  as  shown  in  Figure  1  for  the  colloid  minusil 
in  various  electrolytic  solutions .  The  reason  why  this  curve  is 
of  importance  is  that  it  gives  the  electric  potential  at  the 
surface  of  a  colloidal  particle  from  which  is  calculated  the 
force,  or  more  properly,  the  energy  of  interaction  between  two 
particles  as  they  approach  each  other.  If  the  particles  have  the 
same  surface  potential,  a  curve  representing  the  repulsion  energy 
versus  separation  between  the  particles  is  obtained  as  shown  in 
Figure  2.  When  this  curve  is  combined  with  the  attraction  curve 
due  to  van  der  Waals  or  London  forces,  the  total  interaction 
energy  versus  separation  curve  as  shown  in  Figure  2  is  obtained 
whose  maximum  represents  the  energy  barrier  which  the  particles 
must  overcome  in  order  to  make  lasting  contact.  Since  the  repul¬ 
sion  curve  depends  on  the  value  of  the  surface  potential,  the 
zeta  potential  determines  the  magnitude  of  the  energy  barrier. 
Values  of  the  zeta  potential  around  z^ro  indicate  an  unstable 
system  while  large  positive  or  negative  values  suggest  stability. 

For  the  case  of  the  phosphate  Na^PjO^  examination  of 
Figure  1  shows  that  as  the  concentration  is  increased,  the  zeta 
potential  initially  becomes  more  negative  but  gradually  levels 
off  and  then  becomes  less  negative  until  it  reaches  zero  where¬ 
upon  it  would  no  longer  change  with  concentration.  Similar 
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Figure  1 

CONCENTRATION  OF  ELECTROLYTE 
(from  Riddick,  1968) 


Figure  2 

TOTAL  POTENTIAL  ENERGY  (VT)  DIAGRAM 
FOR  TWO  COLLOIDAL  PARTICLES 
(from  Napper,  1970) 


results  are  obtained  for  KCl,  NaCl,  and  l^SO^.  However,  in  the 
case  of  AlCl^  and  CaC^.  it  is  seen  that  the  zeta  potential 
Initially  becomes  less  negative  upon  Increasing  the  concentration. 
If  the  curves  for  A1C13  and  CaCl2  were  continued  to  higher  con¬ 
centrations,  it  would  be  found  that  the  zeta  potential  would 
cross  zero  and  become  more  and  more  positive  until  it  would  reach 
a  maximum  after  which  it  would  steadily  decrease  to  zero  and  no 
longer  change  with  concentration.  It  is  interesting  to  note  from 
such  data  that  the  same  electrolyte  can  stabilize  or  agglomerate 
a  colloidal  system  depending  on  its  concentration. 

The  second  curious  observation  arises  from  experiments  in 
which  the  adsorption  density  of  ions  adsorbed  onto  a  surface  and 
thus  the  surface  charge  is  measured  for  various  concentrations 
of  an.  indifferent  or  nonadsorbing  electrolyte.  Figure  3  shows 
such  curves  as  determined  for  ferric  oxide  where  the  adsorption 
density  of  hydrogen  and  hydroxide  ions  is  given  as  a  function 
of  pH  and  ionic  strength.  Similarly,  Figure  4  gives  the  surface 
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Figure  3 

ADSORPTION  DENSITY  OF  POTENTIAL  DETERMINING  IONS 
ON  FERRIC  OXIDE  AS  A  FUNCTION  OF  pH  AND  IONIC  STRENGTH: 
TEMPERATURE,  21#C;  INDIFFERENT  ELECTROLYTE  KNO. 
(from  Parks  and  de  Bruyn,  1962) 
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Figure  4 

ADSORPTION  DENSITY  OF  POTENTIAL  DETERMINING  IONS 
ON  SILVER  SULFIDE  AS  A  FUNCTION  OF  THE  pAg 
(from  Freyberger  and  de  Bruyn,  1957) 
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charge  on  silver  sulfide  as  function  of  pAg  and  ionic  strength, 
and  Figure  5  the  charge  on  alumina  as  a  function  of  pH  and  ionic 
strength.  The  observation  to  be  explained  here  is  that  if  the 
surface  is  positively  charged,  an  increase  in  ionic  strength 
makes  the  adsorption  density  or  surface  charge  more  positive. 
However,  if  the  surface  is  negatively  charged,  an  increase  in 
ionic  strength  makes  the  adsorption  density  and  surface  charge 
more  negative.  Also,  the  point  where  the  surface  charge  becones 
zero  is  independent  of  ionic  strength.  Since  the  zeta  potential 
depends  on  both  surface  charge  and  ionic  strength,  it  is  necessary 
to  explain  this  observation  in  conjunction  with  the  previous  one 
in  order  to  obtain  a  good  understanding  of  the  zeta  potential- 
concentration  relationship. 


1 1 T  RESEARCH  INSTITUTE 

6 


Figure  5 

ADSORPTION  DENSITY  OF  POTENTIAL-DETERMINING  IONS 
ON  ALUMINA  AT  25.0°C  AS  A  FUNCTION  OF  pH 
AND  IONIC  STRENGTH,  USING  KCl  AS  SUPPORTING  ELECTROLYTE 

KCI,  NORMALITY:  A - 10"3 ;  O - I0"2;  0 - 10" 1 

(from  Yopps  and  Fuerstenau,  1964) 
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PRESENT  THEORY 


The  present  theory  available  for  analyzing  these  phenomena 
is  the  Gouy-Stern-Grahame  model  of  the  electrical  double  layer 
(Parfitt,  1969).  As  shown  in  Figure  6,  the  double  layer  i*) 
divided  into  three  parts:  (1)  the  region  next  to  the  solid  sur 
face  contains  ions  adsorbed  onto  the  surface,  (2)  the  second 
region  takes  into  account  the  fact  that  counter  ions  which  are 
not  adsorbed  but  are  attracted  to  the  surface  due  to  electro¬ 
statics  cannot  approach  the  surface  closer  than  some  distance 
due  to  their  finite  size,  and  (3)  the  third  part  is  the  diffuse 
layer  as  determined  by  the  Poisson-Boltzmann  equation. 


4»o  4>n  ik, 


Co 

Solid 

Ep 


Diffuse  pa^t 
of  double 
layer 


Figure  6 

MODEL  FOR  THE  ELECTRIC  DOUBLE  LAYER 
AT  THE  SOLID-SOLUTION  INTERFACE 
SHOWING  POTENTIALS  f  AND  CHARGE  DENSITIES  q 
(from  Parfitt,  1969) 
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From  the  condition  that  the  system  be  electrically  neutral 

<!o  +  1b  +  <ld  ■  0  (1) 

where  q  denotes  charge  per  unit  are*  . 

Defining  the  total  capacity  of  the  two  inner  regions  as 

el 

c  ‘  !W3 

where  e  is  the  dielectric  constant  then 


From  the  definition  of  capacity,  one  may  write 


(2) 


and 


♦ 

a 


£ 

0 


+  o. 


(3) 


where  ♦  is  electric  potential.  The  classical  Gouy-Chapman 
equation  is  used  to  relate  q^  and  i|r^ 


_  _  -\|2eRTc  .  ze,|,d 

qd  \Jimttt  Slnh  TF 

where 

R  is  the  universal  gas  constant 
k  is  Boltzmann  constant 
T  is  absolute  temperature 
z  is  valence 
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e  is  charge  of  an  electron 

c  is  molarity  of  electroyte  in  the  bulk  solution 

The  Stern  adsorption  isotherm  is  then  introduced  in  order 
to  relate  the  surface  charge  due  to  adsorption  to  the  bulk  con- 
centration 


zeNQMc  r  ze<|fp  +  $ 

looo  exp  L — e — 


(5) 


where 


Nq  is  the  number  of  adsorption  sites  per  unit  area 


M  is  molecular  weight  of  dispersion  medium 

$  is  the  energy  of  adsorption 

Finally,  the  surface  potential  is  related  to  the  bulk  con¬ 
centration  by  the  thermodynamic  relationship 


♦o 


kT 

e 


In 

Czpc 


(6) 


where  c  is  the  bulk  concentration  for  which  the  surface  charge 
zpc 

is  zero. 


It  is  of  interest  to  list  the  parameters  which  are  not  known 
in  the  Gouy-Stern-Grahame  model.  The  capacities  and  thicknesses 
6f  regions  1  and  2  contribute  the  four  parameters  Cp  0.2*  P  and 
a  while  the  Stern  equation  contributes  the  two  parameters  $  and 
Nq.  One  would  expect  that  with  six  arbitrary  constants,  it 
would  be  possible  to  fit  this  model  easily  to  zeta  potential 
versus  concentration  data  by  judicious  selection  of  the  values 
of  these  constants.  However,  Hunter  and  Wright  (1971)  have 
found  that  this  is  not  necessarily  true  and  had  difficulty  in 
obtaining  agreement  between  the  Gouy-Stern-Grahame  theory  and 
data  when  they  tried  to  apply  it  to  several  sets  of  data  avail¬ 
able  in  the  literature.  In  many  cases,  agreement  could  only  be 
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obtained  when  somewhat  unreasonable  values  of  the  various  param¬ 
eters  were  used,  and  unrealistic  assumptions  were  made.  For 
example,  for  alumina,  it  was  necessary  to  postulate  adsorption 
of  Na  and  Cl  ions  onto  positively  and  negatively  charged  sur¬ 
faces,  respectively.  Such  a  phenomenon  is  known  not  to  occur  to 
any  great  extent  if  at  all.  These  investigators  came  to  the 
conclusion  that  the  Gouy-Stern-Grahame  model  is  inadequate  to 
describe  satisfactorily  the  concentration  dependence  of  the  zeta 
potential  and  that  a  more  sophisticated  model  is  needed,, 

In  the  sections  that  follow,  such  a  model  is  presented.  The 
papers  (see  References)  of  Ottewill  and  his  students,  Parks,  and 
Stumm  and  Morgan  were  of  great  assistance  in  developing  the  con¬ 
cept  of  this  model  and  formed  the  foundations  of  the  theoretical 
work.  It  is  shown  that  this  model  easily  explains  the  experi¬ 
mental  observations  previously  discussed  as  well  as  others  and 
is  in  quantitative  agreement  with  data.  The  reason  why  this 
model  is  more  successful  in  describing  the  concentration  depen¬ 
dence  of  the  zeta  potential  than  was  the  Gouy-Stern-Grahame  model 
is  that  more  physics  was  incorporated  into  the  model.  As  a 
result,  nebulous  quantities  such  as  the  capacities  of  the  inner 
regions  did  not  have  to  be  introduced  in  the  model  so  that  all 
the  parameters  are  of  a  fundamental  and  basic  nature. 
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3. 


SURFACE  CHARGE  DUE  TO  MULT I- COMPONENT  ADSORPTION  FROM 
SOLUTION  (?NTO  A  SOLID  “ 


Surface  charge  on  a  solid  particle  can  be  established  by  the 
adsorption  of  ions  from  solution.  Preferential  adsorption  of 
one  species  over  another  can  result  from  electrostatic  attrac¬ 
tion,  London- van  der  Waals  interaction,  hydrogen  bonding,  and 
chemical  reactions  at  the  surface.  For  example,  cationic  poly¬ 
electrolytes  adsorb  on  negatively  charged  particles  by  electro¬ 
static  attraction  and  non-ionic  and  anionic  ones  by  van  der  Waals 
attraction  and  hydrogen  bonding  despite  the  electrostatic  repul¬ 
sion  in  the  case  of  the  anionic  polymer.  Polynuclear  hydroxo- 
metal  complexes  and  phosphates  adsorb  onto  oxide  surfaces  by 
chemical  surface  reactions.  In  the  sections  that  follow,  these 
adsorption  mechanisms  will  be  discussed,  and  expressions  for  the 
adsorption  isotherms  will  be  derived. 


Under  the  assumption  that  adsorption  of  species  in  solution 
occurs  on  independent  localized  sites  on  the  surfaces  of  a  solid, 
it  has  been  shown  (Fowler  and  Guggenheim,  1965)  that  in  accordance 
with  the  methods  of  statistical  mechanics,  the  chemical  potential 
of  each  adsorbed  species  i  from  solution  onto  a  solid  surface 
may  be  written  as 


where 


(7) 


A  f*  !«* 

is  the  chemical  potential  of  the  i  species 
on  the  solid  at  the  standard  state 
defined  as  half  coverage  of  the  surface 
by  the  ith  species  or  8^  =  0.5 

k  is  Boltzmann  constant 

T  is  absolute  temperature 

8^  is  the  fraction  of  available  sites  occupied 
by  the  ith  species. 
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From  classical  thermodynamics,  the  chemical  potential  of 
the  i*"*1  species  in  the  bulk  solution  is  given  by 

|i1  *  +  kT  In  ,  (8)  1 

t  .  i 

where  .  ' 

is  the  chemical  potential  of  the!  i1"*1  species 

in  solution  at  the  standard  state  defined  1 

as  unit  activity  of  the  it"  species 

x.  is  the  mole  fraction  of  the  ith  species  in 
solution  1  ,  i 

I 

is  its  activity  coefficient. 

At  equilibrium,  the  chemical  potential  bf  the  i4"^  species 
must  be  the  same  in  both  the  solid  and  liquid  phases.  Thus, 
Equations  7  and  8  may  be  equated  to  each  other  to  give  1 


where  -  |i^s  and  is  the  overall  standard  free  energy 

of  adsorption.  This  term  consists  of  the  sum  of  the  specific 
physical  or  chemical  adsorption  energy  and  the  electrostatic  \ 
work  involved  in  transporting  an  ion  frorti  the  bulk  solution  to 
the  charged  surface  and  can  be  written  as  1  1 

MJj0  -  +  Zje*0  ’  , 

I 

where 

t 

!  ,  !  i  | 

§  is  the  energy  per  ion  of  specific  adsorption 
due  to  nonelectrostatic  forces 

is  the  electrical  potential  at  the  surface. 
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The  electrostatic  term  is  obtained  as  follows.  The  work  W 
required  to  transport  a  charge  q  in  an  electric  field  of  strength 
E  from  the  bulk  to  the  surface  is 


i  1  0 

W  -  -  /  Eqdx 


/  q  ||  dX  - 


Since  this  work  is  per  ion,  q  ■  ze  and  W  ■  ze^Q. 

If  is  the  number,  of  molecules  if  the  ic  species  adsorbed 
and  Nj.  the  number  of  adsorption  sites  available  for  it,  then 


I  N?* 

1  1  ■ 

so  that  Equation  9  becomes 


ni  ■  vixiNi' 1  - 1  Sf  bi 


where  for  simplicity  the  term 


'bi  *  exp  “rrj 

has  b^en  introduced.  Equation  11  defines  a  system  of  coupled 
equations  in  the  unknowns  n^ .  Fortunately,  if  the  ratio  of 
Equation  11  for  two  different  values  of  i  is  taken,  n^  can  be 
expressed  in  terms  of  n^  as 

v.x.N.b. 


1  1  M  1 

nj  *  ni 


Then,  substituting  Equation  12  into  11, 

ni  *  .  ViVi  •  ni  1  vjxjbj 

'  1  J 
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or  solving  for 


n. 


vixlNlbl 
1  ^ 


(13) 


ThuSf  in  this  way  it  has  been  possible  to  decouple  Equation  11 
and  solve  explicitly  for  each  n^  in  terms  of  only  v. ,  x. ,  b, , 
and  J  J  J 


The  total  charge  q  on  the  surface  of  the  solid  can  now  be 
determined  by  summing  the  contributions  of  each  species 


<i  ■  I 

i 


(14) 


where 

is  the  valence  of  the  ith  species 

e  is  the  charge  of  an  electron  (1.6  x  10“19  coul) . 

Substituting  Equation  13  into  14,  the  expression  for  the 
charge  becomes 


Zz.eN.  v.x.b. 
7'V 

i  1  +  ^  VjXjbj 
j 


(15) 


Equation  15  is  a  multi-component  version  of  the  Langmuir 
adsorption  isotherm. 


It  is  also  possible  for  the  various  species  to  adsorb  not 
only  on  the  solid  but  also  onto  each  other  so  that  multi-layer 
adsorption  can  occur.  By  analogy  with  Equation  13,  the  number 
of  molecules  of  the  itb  species  which  are  adsorbed  onto  mole- 
cules  of  the  m  species  already  on  the  solid  can  be  written  as 
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(16) 


n  -  VlXl'Vbl- 

iB  1  +  1  Vjbj 
i 

by  considering  nffl  cs  Che  number  of  available  sites  for  formation 
of  a  second  adsorption  layer.  The  notation  for  the  subscript 
MioM  refers  to  the  iC^  species  adsorbed  onto  the  mc^  species  on 
the  solid.  The  additional  charge  on  the  solid  due  to  this  second 
adsorption  layer  is 

^  *  Z  l  Wk»  <17> 

m  k 


or  substituting  Equation  16  into  17 

lq  .  V  y  ,kn.e  VknBbkm 

m  k  1  +I  VJxJb]m 

J 

Then  substituting  Equation  13  into  18 


(18) 


4,  .  y  y  *ta«vkV,km  .  vmxmliapn 

n  k  1  VJXJbJm  1  +  1  vjxJbJ 


v_x_N_b_ 


or 


bq  - 


T7T77V  l  +  \  v  x  b  ^  2l°nVkXkbkm 

1  +  1  JXJbJ  n  1  +  L  vJxJbJn  k 


(19) 


For  the  case  where  each  species  adsorbs  only  onto  species 
of  the  same  kind,  then 


if  k  +  m 
if  k  ■  m 
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Thus,  Equation  19  becomes 


Aq 


1 

1  +  I  vJxJbJ 

j 


I 


2«vxNb 

10  D  D  D  D 


vmxmb- 


1  +  v  x 


HE® 


(20) 


If  the  charges  due  to  the  first  and  second  adsorption  layers 
as  given  by  Equations  15  and  20  are  combined,  the  expression  for 
the  total  charge  on  the  solid  is 


q  j 


1  + 1  vjxjbj 

i 


l 

i 


*ieNivixibi 


1  + 


Vib 


TT 


Vi 


li 


(21) 


In  the  case  where  each  species  has  access  to  the  same  avail* 
able  sites  so  that  -  Hq,  then  Equation  21  becomes 


qT 


V 

1  + 1  Yj' 

j 


-  I  iivtxibi 
’j  i 


,  vixibli 

1  +  tt“v^cf; 


iTii 


(22) 


If  Nq  is  expressed  as  number  of  sites  per  unit  area,  the  the 
units  of  the  charge  is  coulombs  per  unit  area. 
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4.  THE^REL^TIONSHtP  BETWEEN  SURFACE  CHARGE  AND  THE 

4.1  The  Surface  Charge  -  Surface  Potential  Relationship 

The  relationship  between  the  volisae  charge  density  p  and 
the  electric  potential  t  is  governed  by  the  Poisson  equation 

V(«V*)  "  "  4np  (23) 

subject  to  the  boundary  conditions  that  at  the  surface  of  the 
particle  ♦  •  Iq  and  in  the  bulk  solution  f  »  V  ♦  *  0  where  e  is 
the  dielectric  constant.  It  has  been  found  that  the  dielectric 
constant  decreases  with  increasing  electric  field  strength 
according  to  the  empirically  determined  equation  (Booth,  1951) 

*  -  «0  t1  -  *)2  ]  (24) 


where  Cq  is  the  dielectric  constant  at  aero  electric  field 
strength  and 

1x10“ 14  <  b  <  i2xio“  *4  -saL 

voit^ 


Since  large  electric  field  strengths  can  be  realized  in  the 
vicinity  of  a  charged  surface,  it  is  necessary  to  determine 
whether  the  variation  of  the  dielectric  constant  with  field 
strength  is  significant.  If  an  order  of  magnitude  estimate  of 
V ♦  is  considered  as  the  change  in  ♦  across  the  thicl^ess  of  the 
electrical  double  layer,  then  for  the  one»dlmensional  case 

v*  -  &  ~  #  ~  & 

where 

♦g  is  the  surface  p  tentlal 
x  is  distance 


4tt€^ 


I  niori2 

i 
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n.Q  is  the  bulk  number  concentration  of  the  i^ 
species  in  solution. 

The  fact  that  1/*  is  an  estimate  of  the  double  layer  thickness  is 
the  classical  result  of  the  Debye-Huckel  theory  for  electrolytes. 


From  experimental  potential  versus  concentration  curves  (see 
Figures  12  to  17),  maximum  values  of  the  surface  potential  are 
about  150  mv  and  occur  in  the  neighborhood  of  lO"4  moles / l.  Thus 
typical  values  of  the  ionic  strength  and  which  will  give  a 
maximum  effect  are  taken  as  10“4  moles/j l  for  which  1/K  *  300  A 
and  ¥q  ■  150  mv. 


v  ^  150xl0‘\ 

300xl0~^cm 


-  0.50x10 


-5  volt 
cm 


Then,  according  to  Equation  24 

2 

7-  -  1  -  12xl0"14(0.50xl05) 

0 


1  -  3.0xl0‘4 


Thus,  it  is  seen  that  the  effect  of  field  strength  of  the  dielec¬ 
tric  constant  can  be  neglected  so  that  the  dielectric  constant 
can  be  taken  as  being  constant. 

In  order  to  determine  the  charge  density  in  terms  of  the 
potential,  the  Boltzmann  distribution  is  assumed  to  be  applicable 
This  gives  the  number  n^  of  molecules  of  the  ifc^  species  at  a 
position  away  from  the  surface  where  the  potential  takes  on  a 
given  value  as  (Bikerman,  1942;  Schlogl,  1954) 


where  flpj  is  the  volume  of  the  ith  species.  Equation  25  is  a 
more  general  form  of  the  classical  Boltzmann  equation  in  that 
allowance  has  been  made  for  the  space  occupied  by  the  molecules. 
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If  the  ratio  of  Equation  25  is  taken  for  any  two  species, 


then 


n. 


n 


-zie*/kT 
ni0  e  1 

nj0 


(26) 


Therefore,  each  n^  in  Equation  25  can  be  eliminated  to  yield: 

-z^/kT 

- -  (27) 


niOe 


i  r1 

1  +  l 
j 


j0£ 


^ef/FT 


"-Tv 


j 


j  jO 


The  classical  Boltzmann  equation  assumes  that  the  molecules  are 
point  charges  with  zero  volume  so  that 

z^e^ 


nj 


ni0  exp 


i  ^  , 

L  “FT  J 


(28) 


However,  if  the  potential  were  sufficiently  large,  unreasonable 
values  for  the  concentration  near  the  surface  would  be  obtained 
because  there  is  no  limit  on  how  many  molecules  can  be  contained 
in  a  given  volume.  In  order  to  determine  whether  Equation  28  is 
sufficient  or  that  Equation  27  must  be  used,  a  typical  calcula¬ 
tion  of  the  correction  factors  in  Equation  27  will  be  made  for 
a  large  ionic  radius  of  3  A  (1.1  for  N  ,  1.81  for  Cl  ) ,  a  poten- 
tial  of  150  mv,  z  ■  -  1,  and  a  bulk  concentration  of  0.0001  mole/£. 

CP  -  $  t,  (3xl0-8cm)3  -  1.13X10-22 


n«  -  0.0001  6<02xl023  S°lecules  - 1 

U  A  mole 


1000  cm“ 


-  6.02X1016  molec?les 

cnr 
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cpnQ  ■  6.8x10" 

e-ze*/kT  .  150/25  .  403,43 

W0e-ze,|'/kT  -  2. 74x10" 3 

Therefore,  It  Is  seen  that  It  Is  justified  to  use  the  classical 
form  of  the  Boltzmann  equation  as  given  by  Equation  28. 

The  charge  density  can  now  be  determined  by  summing  the 
charges  of  the  molecules 

p  "  I  2ieni  (29) 

and  using  the  Boltzmann  distribution  to  give 


V  ,  ziev  "I 

p  *  1  zienl0  exP  L"  TT  J  <30> 

i 

Substituting  Equation  30  into  23  and  making  use  of  the  fact  that 
the  dielectric  constant  may  be  considered  as  constant,  the 
classical  Poisson-Boltzmann  equation  is  obtained 


v2t 


I  zin10  exP  L"  “PT  J 


(31) 


4.1.1  Planar  Surfaces 

2 

Only  for  the  planar  one-dimensional  case  where  Y7  ^i|i  ■  fL_|. 

dxZ 

can  Equation  31  be  integrated  analytically  (Parsons,  1954).  Such 
a  simplification  infers  that  the  double  layer  thickness  is  small 
in  comparison  to  the  radius  of  the  particle  as  will  be  shown 
later.  The  result  of  a  first  integration  is 


.  +  BnbT 

—  d 


I  ni°  |exp(""ir)  - 
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where  use  has  been  made  of  the  boundary  condition  i|i  ■  ■  0 

as  x  -*  An  expression  for  the  surface*  charge  q  per  unit  area 
in  terms  of  the  surface  potential  can  now  be  obtained  in  the 
following  way.  In  order  to  preserve  electroneutrality,  the  sur¬ 
face  charge  must  be  equal  to  the  negative  of  the  total  charge  ir 
the  solution  so  that 


co 


(33) 


Eliminating  p  in  Equation  33  by  means  of  the  Poisson-Boltzmann 
equation  then  gives 


dx 


X*  OO 

4tt  df 
e  clx 

x-0 


Thus,  substituting  Equation  32  into  34,  the  surface  charge  in 
terms  of  surface  potential  at  x  >  0  is  obtained. 


kTe  V 

Tt r  L  nio  i 

1 exp  1*  “ET  j  *  1 1 

► 

i 

i  j 

4.1.2  Spherical  Surfaces 


(35) 


The  more  difficult  case  of  when  the  assumption  of  a  planar 
interface  is  not  made  will  now  be  discovered.  In  the  case  of 
spherical  particlej,  the  Poisson-Boltzman  equation  as  given  by 
Equation  31  must  be  written  in  spherical  coordinates  as 


*What  is  meant  by  the  surface  from  hereon  is  the  plane  which  the 
charge  resides.  For  example,  in  the  case  of  adsorbed  layers 
the  surface  would  be  considered  as  the  plane  of  these  layers  and 
not  the  surface  of  the  solid. 
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(36) 


^(t2&)  - 
subject  to  the  boundary  conditions 
♦  ■  <|»q  at  r  ■  a 

^  “  cTr  "  0  as  r-oo 

where 

r  is  the  radial  distance  from  the  center  of 
the  particle 

a  is  its  radius . 

If  a  different  radial  distance  is  introduced  which  is 
measured  from  the  surface  of  the  particle  and  scaled  with  respect 
to  the  double  layer  thickness  according  to 

s  -  k  (r  -  a) 

or 

r  "  ^  +  a 


then  Equation  36  becomes 


[(s/lt  +  a)2  37*7*7] 


fr'p(t) 

e 


or 


(s 

subject  to 


(37) 


♦  ■  at  s  ■  0 


♦  ■  ■g-j-  "0  at  s  00 
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Since  '1 1  varies  mainly  in  the  double  layer  or  between  r  ™  a 
and  r  -  a  +  1/K,  the  condition  at  infinity  can  be  replaced  by  a 
finite  distance.  For  example,  numerical  calculations  showed  that 
a  distance  of  three  double-layer  thicknesses  represented  infinity 
very  well 


* 


df 

3s 


0 


at 


s  —  3 


The  condition  for  which  the  solution  for  the  planar  inter- 
face  approximation  is  valid  can  be  easily  obtained  from  this 
analysis.  Equation  37  shows  that  as  Ka  becomes  large,  it  will 
dominate  in  the  term  "s  +  Ka."  Since  s  is  on  the  order  of  unity, 
it  is  seen  that  for  values  of  Ka  much  greater  than  one,  Equation  37 
reduces  to  that  for  the  planar  case 


In  order  to  calculate  the  surface  charge,  one  proceeds  as  in 
the  planar  case.  The  total  charge  Q  on  the  particle  is 


oc 


oo 


Q  =  -  /  4nr2pdr  •=  /  “jr  (*2  $jr ) 


a 


dr 


oo 


f  d 

J  IF 


M) 


oo 


dr 


cr2  31 

er  3F 


a 


r=a 


Thus,  the  charge  q  per  unit  surface  is 


Q/4TTa' 


e  df. 
vlr 


r=a 


(38) 
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For  the  spherical  interface,  an  exact  analytical  solution 
is  not  possible.  However,  if  the  charge  density  p  is  expanded 
in  series  about  t  ■  0  and  only  the  linear  term  retained,  a  solu¬ 
tion  valid  for  ze^/kT  <  1  or  <  25  mv  can  be  obtained  as  given 


by  the  Debye-Huckel  theory  (Tanford,  1961). 
linearized  Poisson-Boltzmann  equation 

The  solution  of  the 

1  d 
7Hr 

(<24f)  ■ 

K2)| i 

(39) 

subject  to 

*  = 

*0  at 

r  ■  a 

*  - 

O 

a 

■Si'S 

as  k 

yields 

*  “ 

*0  £  e*P  L- 

K  (r  -  a)  ] 

(40) 

q  - 

ei(i0(l  +  Ka)/4na 

(41) 

For  practical  purposes,  this  solution  is  not  very  useful  because 
the  magnitude  of  the  surface  potential  required  for  stabilization 
of  colloidal  systems  greatly  exceeds  25  mv  by  factors  of  2  to  6. 
Therefore,  the  case  of  a  spherical  interface  must  be  solved 
numerically  by  computer. 

Mathematical  problems  of  this  type  known  as  two  point  bound¬ 
ary  value  problems  because  the  conditions  on  are  known  at  the 
particle  surface  and  in  the  bulk  solution.  The  problem  would  be 
relatively  simp  3  if  both  the  potential  and  its  derivative  were 
known  at  the  surface  because  then  the  non-linear  differential 
equation  could  be  integrated  numerically  like  an  initial  value 
problem. 
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However,  since  the  problem  is  to  calculate  the  surface 
charge  due  to  a  given  surface  potential,  such  a  procedure  cannot 
be  followed  because  the  first  derivative  is  proportional  to  the 
surface  charge  and  thus  not  known.  Fortunately  on  another 
project,  the  author  has  already  developed  a  computer  program 
which  solves  the  more  general  multi-point  boundary  value  problem. 
The  basis  of  the  numerical  scheme  that  is  employed  by  this  pro¬ 
gram  is  as  follows.  In  order  to- take  advantage  of  the  simplicity 
of  solving  an  initial  value  problem,  an  initial  guess  of  the  first 
derivative  or  surface  charge  is  made.  This  allows  the  differ¬ 
ential  equation  to  be  integrated  numerically  from  the  surface  to 
the  bulk.  At  this  point,  the  program  checks  to  see  whether  the 
calculated  conditions  for  the  bulk  solution  are  in  agreement  with 
those  specified.  If  they  are,  the  initial  guess  was  correct,  and 
the  problem  is  solved.  However,  if  they  are  not,  the  program 
calculates  a  first-order  correction  to  the  initial  guess  by  means 
of  a  combination  of  the  Newton-Raphson  method  and  a  gradient 
search  method  which  insures  convergence.  This  corrected  value 
of  the  first  derivative  is  then  used  to  integrate  the  differ¬ 
ential  equation  a  second  time.  Again  the  bulk  conditions  are 
checked  to  see  if  they  are  satisfied  yet.  This  interactive  pro¬ 
cedure  continues  until  the  bulk  conditions  are  satisfied. 

A  comparison  is  given  in  Table  1  of  the  surface  charge  as 
calculated  for  a  planar  and  spherical  interface  where  the  surface 
potential  and  ionic  strength  are  the  same.  It  is  seen  that  the 
charge  for  the  planar  surface  is  smaller  than  that  for  the 
spherical  surface  but  that  as  Ka  increases  the  difference  between 
them  gradually  decreases  as  expected. 
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Table  1 

COMPARISON  OF  SURFACE  CHARGE  FOR  PLANAR 
AND  SPHERICAL5  SURFACES* 


Ka 

*0 

(mv) 

^planar 

2 

(coul/cmf) 

| 

^sphere 

2 

.(coul/cm  V 
-5 . 163xlf0“7 

0.660 

-80 

-2 . fSOxlO"7 

0.662 

-80 

-2.667xlO-7 

•  ! 

-5.l70xl0‘7 

0.664 

-68 

-2.061xl0"7 

-4.261x10"? 

0.670 

-40 

-1.024x10" 7 

-2 .386xl0-7 

0.682 

-16 

-3.827xl0"8 

-9.376xl0"8 

0.701 

12 

2.933xl0"8 

71  101xl0"8i 

0.711 

20 

5.039xl0"8  ' 

1. 196xl0"7 

| 

0.757 

48 

1.448xl0"7 

3  .,055x10"  7 

0.885 

72 

2.996x10-7 

5.282xlQ"7 

0.990 

88 

4 . 726xl0"7 

i 

7. 378x10" 7 

*Data  taken  from  the  zeta  potential  versus  concentration  curvfe 
for  DTAB. 

a  ■  200  A  (see  Figure  14) .  1 


4.2  The  Zeta  Potential 

-  ,  , 

In  electrophoresis  experiments,  the  mobility  of  a  charged 
particle  is  measured  by  observing  the  velocity  at  which  it  moves 
in  an  imposed  electric  field.  This  mobility  catf  then  be  used  to 
calculate  the  surface  potential  as  discussed  in  detail  in  the 
literature  review.  Although  general  equations  were  presented;, 
they  were  all  based  on  the  assumption  that  the  viscosity  in  the 
vicinity  of  the  particle  is  equal  to  that  in  the  bulk  solution. 
However,  since  viscosity  varies  with  electric  field  strength, 

5  1  ■ 

i  « 

i 
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this  assumption  may  not  be  valid.  Following  Lyklema  and  Over- 
beek  (1961),'  this  problem  will  be  investigated  by  modifying  the 

f  !  • 

Smoluchbwski  equation. 

i 

i  If  the  thickness  of  the  electrical  double  thickness  is  much 
smaller  than  the!  radius  of  the  particle,  the  liquid-solid  inter- 
face  can  be  treated  as  a  planar  surface.  Then  a  volume  element 
in  the  solution  near  the  moving  charged  surface  of  thickness  dx 
and  unit,  area  (see  Figure  7)  experiences  forces  due  to  the 
imposed  electric,  field  ahd  viscosity.. 


i  i 


Ax 


f 


EPdx 


dv 

u  35 


777777/777777 

'/  II 1 1  hi:  i!l !  h 


in 1 


U  <- 


->  Surface 


Figure  7 

THE  FORCES  ACTING  ON  A  VOLUME  ELEMENT 
NEAR  A  CHARGED  SURFACE 


At  steady  sta^e,  the  net  sum  of  the  forces  must  equal  zero 
so  that 


*.  EpAx  + 


q  dv 
r  ?  lx 


x+Ax 


(42) 


where 


1 

E  is  the  electric  field  strength 

1 

i  1 
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|i  Is  viscosity 
v  is  velocity 
x  is  distance. 


Taking  the  limit  Ax  goes  to  zero,  Equation  42  gives 


Eliminating  p  by  means  of  the  Poisson  Equation  43  becomes 


(43) 


d 

dv  \ 

E 

d 

/.  d*\ 

as 

(M«) 

as 

('=) 

(44) 


Equation  44  is  integrated  from  the  plane  of  shear  where  the 
velocity  is  equal  to  that  of  the  particle  u  to  the  bulk  solution 
where  the  velocity  and  potential  gradients  are  zero.  Due  to  the 
fact  that  the  viscosity  increases  with  increasing  field  strength, 
it  is  possible  that  the  shear  plane  will  not  coincide  with  the 
solid  surface.  A  first  integration  of  Equation  44  gives 


dv  m  sE  dv 
<Ix  4nu  <Tx 


and  a  second  one 


u 


0 


Since  mobility  U  is  defined  as 


U 


u 

I 


Equation  45  becomes 


U 


i 


(45) 


(46) 


Experimental  values  mobilities  are  generally  reported  in 
terms  of  the  zeta  potential  C  which  is  the  electric  potential  at 
the  shear  plane.  It  is  assumed  that  the  dielectric  constant  and 
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viscosity  remain  constant  at  their  bulk  values  from  the  bulk  to 
the  shear  plane  and  then  from  this  point  to  the  surface  the 
viscosity  becomes  very  large  such  that  a  viscous  immobile  layer 
is  formed. 


U 


By  writing  Equation  46  in  the  form 


/ 


shear 

plane 


oo 


dx+  f 

shear 

plane 


(47) 


it  is  seen  that  under  these  assumptions  the  second  integral 
vanishes  and  e  and  n  are  equal  to  their  bulk  values  Cq  and  Pq  in 
the  first  integral  so  that 


shear 


(48) 


If  Equations  46  and  48  are  combined,  then  tg  can  be  related  to 
the  measured  C  according  to 


C 


(49) 


As  noted  earlier,  the  dependence  of  e  on  field  strength  is  given 
by 

c  -  c0  [l  -  B  ]  (50) 

where  1x10'*^  <  B  <  12x10*^  cm  n 
~  ~  volt 
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It  has  also  been  found  that  the  dependence  of  viscosity  on 
field  strength  can  be  expressed  by  the  empirical  equation 
(Hunter,  1966) 

w  ■  “0  [l  +  £  ({&)  2]  (51) 

2 

where  1x10" 13  <  f  <  2x10“ 13  — ■ 

voltz 


Substitution  of  Equations  50  and  51  into  49  yields 


(52) 


An 

d|jf 

Hx 

of 


order  of  magnitude  estimate  of  ^jj-  was  made  earlier  which  gave 


~  0.50xl03  •  Inserting  this  result  and  the  maximum  value 

f  into  Equation  52  gives 


C  -  /  °  1---1-^.1.(>.'14(0-50x105^  d,  .  f  *°  1  -  3.0x10-*  df  ..  t 

0  1  +  2x10"13(0.50x105)  0  1  +  5*0xl° 

This  result  indicates  that  the  zeta  surface  potentials  may 
be  for  ail  purposes  considered  equivalent.  Numerical  integrations 
of  Equation  52  where  ^  was  calculated  from  the  solution  of  the 
Poi88on-Boltzmann  with  an  electric  field  dependent  dielectric 
constant  were  made  by  Hunter  (1966)  whose  results  are  in  agree¬ 
ment  with  the  above  conclusion. 
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5.  EXPLANATION  OF  ZETA  POTENTIAL  VERSUS  CONCENTRATION  CURVES 


It  is  now  possible  to  explain  the  shape  of  zeta  potential 
versus  concentration  curves.  From  the  mathematical  theory  pre¬ 
sented,  a  pictorial  description  of  the  electrical  double  layer 
can  be  constructed  as  shown  in  Figure  8.  A  particle  ionersed  in 
aqueous  solution  can  develop  a  surface  charge  by  adsorbing  ions 
denoted  as  potential  determining  ions  onto  its  surface.  As  a 
result  of  this  surface  charge,  ions  of  the  opposite  charge  will 
be  attracted  to  the  surface  while  those  of  the  same  sign  will  be 
repelled.  In  this  way,  a  concentrated  layer  of  ions  oppositely 
charged  to  the  surface  is  formed  around  the  particle.  Assuming 
that  these  ions  do  not  adsorb  onto  the  surface,  one  denotes  them 
as  counter  ions.  The  thickness  of  this  double  layer  is  such 
that  this  layer  contains  a  sufficient  number  of  counter  ions  to 
neutralize  the  surface  charge  because  the  system  must  be  elec¬ 
trically  neutral.  Thus,  as  the  bulk  concentration  of  counter 
ions  increases  with  that  of  the  potential  determining  ions 
remaining  constant,  this  thickness  decreases  because  less  volume 
is  required  to  contain  enough  counter  ions  to  neutralize  the 
surface  charge.  From  an  order  of  magnitude  analysis  of  the 
Poisson  equation 


Lii 

(Ax)2 


P 


Ax 


* 


0 


/v> 


Ax 


Therefore,  the  surface  potential  or  the  zeta  potential  follows 
the  behavior  of  the  double  layer  thickness  so  that  they  both 
approach  zero  together. 

From  these  ideas,  the  concentration  dependence  of  the  zeta 
potential  can  be  easily  explained.  With  reference  to  Figure  1 
and  the  Na2P20^  curve,  when  the  concentration  is  initially 
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Figure  8 

CONCEPT  OF  THE  ZETA  POTENTIAL 
(from  Riddick,  1968) 


Layer  of 
Adsorbed 
Ions 


Extent  of 
Diffuse  Layer 
of  Counterions 


« — 


Electric 
Potential 
Surrounding 
the  Particle 


Plane  of 
Shear 

Bulk  of 
Solution 


Concentration 
of  Negative  Ions 


Concentration 
of  Positive  Ions 


1 1 T  RESEARCH  INSTITUTE' 
11 


m 


increased,  the  potential  determining  ions  P207"^  are  adsorbed 
which  decreases  the  magnitude  of  the  surface  charge  and  zeta 
potential  while  the  increase  in  ionic  strength  or  the  concentra¬ 
tion  of  counter  ions  Na"*"  decreases  it.  However,  at  this  stage, 
the  adsorption  effect  is  dominate  so  that  the  zeta  potential 
decreases.  However,  eventually  all  of  the  surface  sites  become 
occupied  so  that  then  the  ionic  strength  effect  prevails,  and  the 
zeta  potential  increases  to  zero.  In  the  case  of  AlCLj,  the  same 
reasoning  can  be  applied  ^ith  the  modification  that  positive  ions 
are  being  adsorbed.  As  shown  later,  Al+3  ions  are  not  adsorbed 
but  rather  its  hydroxo  complexes. 


* 
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6. 


SURFACE  CHARGE  DUE  TO  CHEMICAL  REACTIONS  AT  THE  SOLID 
SURFACE - 


The  surface  charge  on  a  solid  can  be  acquired  as  a  result 
of  chemical  reactions  at  the  surface.  Such  reactions  can  occur 
because  solid  surfaces  contain  ionizable  functional  groups  such 
as  -OH,  -COOH,  -OPO^^.  Thus,  the  surface  charge  depends  on  the 
extent  of  ionization  of  these  groups  which  in  turn  depends  on  the 
pH  of  the  solution.  For  example,  in  an  aqueous  environment,  a 
hydroxylated  surface  is  formed  on  oxidic  solids  due  to  hydrolysis 
of  the  oxide.  The  surface  charge  then  develops  by  means  of  dis¬ 
sociation  of  the  surface  hydroxide  groups  according  to  the 
reactions  (Parks,  1967) 


-M-OH(S)  ~ 

-M-0“(S)  +  H+(aq) 

(Rl) 

-M-OH(S) 

-M+(S)  +  OH“(aq) 

(R2) 

\ 

h20 

M0H2+(S) 

where  M  denotes  the  metal  and  S  the  surface.  Such  a  mechanism  is 
depicted  in  Figure  9  for  the  origin  of  the  electrical  charge  at  ' 
the  quartz  surface  in  aqueous  solutions.  This  reaction  scheme 
demonstrates  the  amphoteric  behavior  of  the  surface  in  that  it 
can  act  either  as  an  acid  by  giving  up  protons  (Rl)  or  as  a  base 
by  accepting  protons  (R2) .  At  low  pH  values,  the  backward  reac¬ 
tion  of  Rl  and  the  forward  reaction  R2  will  be  favored  so  that 
the  surface  charge  will  be  positive;  on  the  other  hand,  at  high 
pH  values,  the  reverse  will  be  tr’ie  so  that  a  negatively  charged 
surface  results.  At  some  intermediate  pH  value  denoted  as  the 
isoelectric  point  or  zero  point  of  charge,  the  surface  charge 
will  be  zero. 

In  addition  to  charging  due  to  dissociation  of  surface 
groups,  charged  hydroxo  complexes  which  are  formed  as  a  result 
of  the  solubility  of  the  solid  can  be  adsorbed  onto  the  solid 
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Figure  9 

MECHANISM  FOR  THE  ORIGIN  OF  THE  ELECTRICAL  CHARGE 
AT  THE  QUARTZ  SURFACE  IN  AQUEOUS  SOLUTIONS 
(from  Gaudin  and  Fuerstenau,  1955) 

surface.  A  reaction  scheme  for  this  phenomenon  can  be  written 
as  (Parks,  1967) 

M203(S)  +  H20^  2M(OH) 3 (aq)  2M(0H)^"r‘(aq)  +  2  (3-n)OH"  (aq) 

2M(OH)p  n(S)  (R3) 

Possible  reasons  why  the  metal  hydroxo  complexes  are  strongly 
adsorbed  at  the  surface  have  been  given  by  Stumm  and  O'Melia 
(1968):  (1)  hydrolyzed  species  are  larger  and  less  hydrated  than 
non-hydrolized  species;  (2)  the  replacement  of  aquo  groups  by 
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hydroxo  groups  in  the  coordination  sheath  of  a  metal  ion  may 
render  a  complex  more  hydrophobic  by  reducing  the  interaction 
between  the  central  metal  ion  and  the  remaining  aquo  groups. 

This  reduction  in  solvent-hydroxo  complex  interaction  can  then 
enhance  the  formation  of  covalent  bonds  between  the  metal  atom 
and  specific  sites  on  the  solid  surface  by  reducing  the  energy 
necessary  to  displace  water  molecules  from  the  coordination 
sheath.  Also  for  polyhydroxo-polymetal  species,  adsorption  is 
especially  enhanced  because  more  than  one  hydroxide  group  per 
molecule  can  become  attached  at  the  interface.  In  addition  to 
metal  ion  hydrolysis  products,  most  species  which  contain  hydrox“ 
ide  groups  in  their  ionic  structure  including  both  c  'tionic  and 
anionic  hydroxo  complexes  adsorb  at  solid-liquid  interfaces. 

For  example,  in  tlfe  ca3e  of  phosphates,  ligand  exchange  can  occur 
according  to  (Stumm  and  Morgan,  1970) 

-M-OH(S)  +  HP0^2(aq)^=^  -M-0P03H"(S)  +  0H“(aq)  (R4) 

Hydrous  metal  oxides  or  hydroxides  can  also  interact  with 
cations  as  a  coordination  with  electron  acceptors  as  given  by 
(Stumm  and  Morgan,  1970) 

-M-OH(S)  +  A+Z(aq)^=*  -M-0A+Z(S)  +  H+(aq)  (R5) 

It  is  difficult  to  say  whether  these  reactions  are  a  result 
of  coordination  by  covalent  bond  formation  or  electrostatic  ion 
exchange.  However,  some  indication  can  be  obtained  from  the 
electronegativity  by  using  the  following  empirical  equation  to 
calculate  the  ionic  character  of  the  bond  A-B  (Moore,  1964) 

%  ionic  character  -  16  |  XA  -  Xfi  |  +  3.5(XA  -  Xfi)2 

where  X  is  the  electronegativity  given  in  Table  2.  For  example, 
in  the  case  of  the  Al-0  bond 

%  ionic  character  «  16  | 3.5  -  1 .47 |  +  3. 5(3. 5  -  1.47)2  «  46 
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The  covalent  character  of  this  bond  gives  some  insight  into  why 
it  is  observed  that  ligands  containing  oxygen  such  as  phosphates 
tend  to  adsorb  onto  metal  oxides. 

As  a  rule  of  thumb,  surface  coordination  by  covalent  bonding 
is  favored  by  small  differences  between  the  electronegativities 
of  the  two  bonding  species  while  large  differences  indicate  an 
electrostatic  interaction. 

Other  examples  where  solutes  become  bound  to  solid  surfaces 
are  (Stumm  and  Morgan,  1970) 

FeOOH(S)  +  H+(aq)  Fe(0H)£(S) 

FeOOH(S)  +  0H~(aq)  +  H20  ^  Fe(0H)^(S) 

V 

FeP04(S)  +  H2P0“(aq)^=^  Fe(HP04)"(S) 

Si02(S)  +  0H“^=^  Si02(0H)"(S) 

AgBr(S)  +  Br“^=^  AgBr2(S) 

Cu(S)  +  2H2S^=±  Cu(SH)“2(S)  +  2H+(aq) 

S(S)  +  S~2  (aq)  S^CS) 

j-9  x  r  ,-n+2m 

R(C00H)n(S)  +  mCa  Z  ^=^R-[(C00)nCamJ  (S) 

+  nH+  (aq) 

Mn02*H20(S)  +  Zn+2  Mn00H0Zn+(S)  +  H+(aq) 
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It  is  of  interest  to  note  that  when  adsorption  results  due 
to  surface  chemical  reactions,  it  can  occur  against  electro¬ 
static  repulsion  forces  because  the  energy  released  due  to  the 
chemical  reaction  can  outweigh  the  electrostatic  work  involved 
in  bringing  the  adsorbant  from  the  bulk  to  the  surface. 

The  composition  and  properties  or  simply  the  reactivity  of 
the  solid  plays  an  important  role  in  determining  the  surface 
reactions  which  occur  and  to  what  extent.  For  example,  Table  3 
shows  that  the  isoelectric  point  for  aluminum  oxide  can  vary  from 
pH  5  to  pH  9.25  as  a  result  of  different  modifications  of  the 
material  and  different  active  forms  of  the  same  modification. 

A  summary  of  the  transformations  of  the  various  crystalline 
modification  of  aluminum  oxide  and  hydroxide  is  shown  in 
Figure  10. 


Table  3 

ISOELECTRIC  POINTS  OF  ALUMINUM  OXIDE 
(from  Parks,  1965) 


Material 

zpc  (pH) 

a-Al203 

5.0,  6.6,  9.2 

V-AI2O3 

8.0 

a-AlOOH 

7.7,  9.1 

a-Al(0H)3 

5.0 

v-Al  (OH)  3 

9.25 
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Gibbsite  (hydrargillitc) 
M(OH),  (table 


Dehydrate 
above  1511* 


Rr  hydrate 
above  140* 


X-Al,0,  — ►  x-AI,(J3 


Rrliyrirair 
hrlnw  l|»“ 


.  Btihmite 
AIO(OH) , 


Wairr  under  pirnur. 


(71° 


Reliyilr»lf^^>  450“, 

J-AI.O, 


V 

(J-Al.O, 


Diaipore 

AlO(Ori) 

|  4i''' 
«-Ai,d, 
(table 

/ 


0-AI,O, 


Itaytrile 
AI(C)II), 
me  tasiable 


Figure  10  ;  n 

SUMMARY  OF  THE  TRANSFORMATIONS  OF  THE  VARIOUS 
CRYSTALLINE  MODIFICATIONS  OF  THE 
ALUMINUM  OXIDES  AND  HYDROXIDES 
(from  Reriy,  1956) 


Active  forms  of  a  material  are  very  fine  crystalline  preci¬ 
pitates  with  a  disordered  lattice  and  are  generally  formted 
incipiently  from  strongly  oversaturated  solutions.  Active  forms 
convert  into  more  stable  or  inactive  f:orms  although  such  trans¬ 
formations  may  proceed  slowly.  Inactive  or  a£ed  forms  have 
ordered  crystals  and  are  thus  less  soluble  than  active  forms.  ■ 
Hydroxides  and  sulfides  are  examples  of  materials  which  occur 
in  amorphous  and  several  crystalline  modifications.  In  the 
case  of  iron  hydroxide,  the  amorphous  active  form1  can  age  to  the 

inactive  form  and  to  more  stable  modifications,  according  to  the 
scheme.  ’  i 
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(amorph>Fe (OH) 3  (active) 


(amorph)FeOn/2(OH)3_n  (inactive) 
a-FeOOH  (active) 


ot.-Fe203 


6*1  Surface  Charges  Due  to  Adsorption  of  Complex  Ions 
Derived  from  Electrolytes 

It  has  been  shown  in  many  studies  that  multivalent  metals 
ions  react  with  water  or  hydrolyze  to  produce  a  variety  of  com- 
ions.  For  the  case  of  trivalent  metal  ions  such  Al+3  and 
Fe,  ,  the  successive  acid-base  reactions  are  (Sullivan  and 
Singley,  1968) 

MC’HjO)^3  =t±  M(H2P),OH+2  +  H+ 

I  i  ! 

1  ' 

^=±  M(H20)4(0H)J  +  2H+:^  M(0H)3(H20)3  +  3H+ 

j  . 

7=^  M(0H)4,(H20^2  +  4H+^z±  M(0H)5H20“2+5H+ 


M(OH)“3  +  6H+ 


(R6) 


For  .bivalent  ions, *  such  as  Zn+2  and  Cu+2,  a  similar  reaction 
scheme  is  followed  except  that  only  four  water  molecules  instead 
of  six  can  be  coordinated  with  the  metal  ion. 

i  '^n  addition  to  these  hydroxo  complexes,  polynuclear  hydroxo 
complexes  are  formed  through  polymerization  of  the  hydroxo  com¬ 
plexes  by  condensation  reaction  such  as 


.OH. 

2M(H20)50H+2  ^=^[ch20)4m^'  \m(h20)4] 

X0H'/ 


+4 


+  2H20 


'  ,  (R7) 
This  dimer  can  then  hydrolyze  to  provide  hydroxo  groups  so  that 
anbther  condensation  reaction  can  occur.  For  aluminum,  the 
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accepted  general  form  of  such  polynuclear  species  Is  Aln(OH)^j^x 
(Brossett,  1954).  The  main  complexes  with  this  formula  which 
have  been  suggested  by  various  Investigators  are  Alg(OH)^ 
(Brossett,  1954)  and  Alg(OH)^Q  (Matijevic,  1961)  and  others  such 
as  A1?(0H)^  (Hsu,  1965)  and  Al^OH)^  (Rausch,  1964).  These 
polynuclear  complexes  have  been  proposed  to  fit  into  the  hydroly¬ 
sis  reactions  as  intermediates  between  AlO^O^OH)^  and 
AlfaH)^ (HjO)^  as  shown  In  Figure  11  (Stumm  and  Morgan,  1962). 
Equilibrium  constants  for  hydrolysis  and  complex  formation  for 
aluminum  and  iron  are  given  in  Table  4. 


Figure  11 

STEPWISE  CONVERSION  OF  THE  TRI-POSITIVE 
ALUMINUM  ION  TO  THE  NEGATIVE  ALUMINATE  ION 
(from  Stumm  and  Morgan,  1962) 
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Tabic  4 


HYDROLYSIS  AND  COMPLEX  FORMATION  EQUILIBRIA 
OF  IRON  AND  ALUMINUM 
(from  Stuam  and  Morgan ,  1970) 


N» 


1 

2 
J 

4 

5 

6 
7 
C 
♦ 
10 

11 

12 
IJ 

14 

15 

16 
17 
IS 
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7n 
21 
22 
a 

24 

75 

7b 

27 


uw 

tHWita 


M 
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Al  ‘  *  *  V  IM)  -  [AKIII)'*  4  II* 

2\l***  4  211/1  -  [AMitlll,}*  4  211* 
Al*'*  4  JIliO  -  AI((HII,  M  4  ‘If 
AUOII),  (4  4  11(0  -  lAI(OII|,V  4  II* 
Ml"  *4  Tfill.O  «  r.\U(Otl)n>'  4  20  11* 
Ml***  4  1511,0  -  fAMOII).,)***  4  ISII* 
IV*  4  SO.**  -  FcM i. ' 

If*  f  111*0.—  -  11*111-0.)* 

Kf  *  4  1*0. —  -  1.1*0.  (.) 

Al"  I  1*0, —  -  AllSI,  (.) 
if*  4  iw- -  [rri',0,)- 

Al***  4  I'Oi**  -  tAII*iOi)* 

If  *  H  -CIVIT* 

Al***  4  7  -  -  [All  )** 

IV  **  4  AfriJtr'  -  fIV  .Wuir)*  * 

Al***  4  I  A*vi  Me*  -  f  Al  (  Vruii)  ] 

If  4  JOuilMr-  *  ••  [IV  (Ou.l.t,4.) - 

Al*  *  *  4  J  Ovilatr  -  f  MllluUirt  ]  *  * 
If  rl.**4(*  -  [IV  |1h-.,4m.  )*♦ 
IV*'*  4  IVlrjIe*  •  {lr filMcV 
le**'  •*  S.lioUlc**  -  (IVSil*)lMeV 
Al  **  4  Nili»  >l4»r~  -  [Al  Silw>hu  J* 


4  H.8> 
4  II* 


IM*<  » 


-  217 

-  I  7 
1.46 

-  6 
-155 

-  5.01 

-  6  27 

-  9.1 

-ll.W 

165 

7.21 

•  4 

219 

22 


5.17 
6  IJ 


7.1.9 
16  J 
6.2 
75 

I5J17 


Compounds  Chat  are  commonly  used  to  supply  aluminum  and 
ferric  ions  In  solution  for  the  purpose  of  stabilization  and 
destabilization  of  colloids  by  hydrolyzed  metal  ions  are  FeCl3, 
Fe(C104)3,  Al2(SO^)3,  and  alums  such  as  Rj^O^Rj^CSO^O^O^ 

or  RIRIII(SO^>2(H20)^2  where  R1  denotes  an  univalent  metal  such 
as  sodium  or  potassium  and  R111  is  either  Al+3  or  Fe+^. 

In  order  to  demonstrate  how  these  ideas  can  be  put  into  a 
quantitative  form,  let  us  consider  the  hydrolysis  of  a  trivalent 
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metal  without  polynuclear  complex  formation  which  can  be  repre¬ 
sented  as 


kl 

M(H20)7-i(0H)i-l  F=i,,(H2O)6.i(a0]"i  +  H+  i-1,2 . 6 


species  i 


species  1+1 


The  equilibrium  relationships  are  then  given  by 


i-1,2 . 6 


(R8) 


(53) 


where  c  denotes  concentration. 

If  it  is  assumed  that  the  amount  of  metal  ions  adsorbed  onto 
the  solid  is  small  in  comparison  to  the  amount  initially  supplied 
to  the  solution  through  an  electrolyte,  a  mass  balance  on  M 
neglecting  the  adsorbed  species  yields 


7 

l  cj  *  °0 

j=l 

where  cQ  is  the  initial  or  total  concentration  of  M. 
Since,  from  Equation  53 


ci+l 


ci 


one  may  proceed  by  induction  to  show  that 


(54) 
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c2 


C1 


c3 


i  k 

TT  K, 

_  1-0  J 

ck+i  "  T TTTcl 

.  lH  J 

i 

i 

.  K1K2K3K4K5K6 

7  [h+]6  Cl 

where  Kq  —  1  and  tt  denotes  the  product  sign. 
Substituting  Equation  55  into  54 
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so  that 


u  h  rw. 

y  .i-°  ,J 

A 

j_U 

Thus,  the  concentration  of  the  k+lL  species  is  obtained  by 
substituting  Equation  56  into  55 


ck+l 


k=0, 1 . 6  (57) 


In  order  to  calculate  the  surface  charge  due  to  adsorption 
of  the  hydroxo  complexes,  Equation  15  derived  earlier  for  surface 
charging  due  to  multi-component  single  layer  adsorption  is  used. 
It  has  been  suggested  that  the  specific  adsorption  energies  for 
all  the  complexes  are  approximately  equal  (Parks  and  de  Bruyn, 
1962). 

Complex  formation  is  not  restri'ted  to  hydroxides  but  other 
groups  such  as  phosphates  can  also  undergo  complex  formation. 

When  a  phosphate  such  as  KHgP04  is  put  into  solution,  the  follow¬ 
ing  equilibria  reactions  occur 

ki  kg 

H3P04  —  v  H2P°4‘  +  H+  . . *  HP04"2  +  2H+ 

species  1  species  2  species  3 

-  P04“3  +  3H+  (R9) 

species  4 
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The  concentrations  of  the  various  species  are  given  by  Equation  57 
where  k  -  0, 1,2,3  and  the  summation  in  the  denominator  ranges 
from  0  to  3. 

Equilibrium  constants  for  the  hydrolysis  reactions  of  various 
phosphates  are  given  in  Table  5. 


Table  5 

ACIDITY  AND  HYDROLYSIS  OF  PHOSPHATES  AND  METAL  IONS 
(from  Stumm  and  Morgan,  1970) 


No. 

Equilibrium 

Log  Eiouilibrium 
Constant,  25  C 

S 

H  .HO,  » 

iijPOi "  +  ir 

-2.2 

11*1*0,  - 

=  HP<V  -1  H  - 

-7.0 

10 

IIPO,=  -  : 

-  pov  +11' 

-12.0 

II 

ll.,IU)7 

=  H,p2<v  -  i-  h  *• 

-2.4 

12 

iww 

-  hp2<v  +  h  + 

-6.6 

1.1 

lllW- 

*  Pi<V  +  II  4 

—  9.3 

II 

H.I’Ao2 

•---  H»PsOl0s*  +  H  + 

-2.3 

15 

ll-.PaO,,,' 

-  HP.,O,04-  +  H + 

-6.5 

l<> 

HP,O10' 

-  P:.O10r'  -f  H  * 

-9.2 

17 

HPAr 

=  PO,' 

-2.1 

These  phosphates  can  then  adsorb  onto  hydrous  metal  oxides 
or  hydroxides  by  ligand  exchange  according  to  (Stumm  and  Morgan, 
1970) 

-M-OH(S)  +  HP04"2(aq)  -M-OPO3H' (S)  +  OH"(aq)  (RIO) 

Although  OH  has  a  greater  affinity  for  metal  ions  such  as 
Fe  than  the  phosphates  do,  the  adsorption  of  the  phosphates  is 
favored  by  reduction  in  pH.  However,  this  reduction  is  limited 
by  the  hydrolysis  equilibria  of  the  phosphates  because  as  pH  is 
decreased  the  fraction  of  the  ionized  phosphates  will  decrease. 

In  order  to  calculate  the  surface  charge  due  to  phosphate  adsorp¬ 
tion,  Equation  15  for  multi-component  adsorption  may  again  be 
used. 
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6.2  Surface  Charge  Due  to  Adsorption  of  Complex  Ions 

Derived  from  the  Solid 

Due  to  the  solubility  of  the  solid,  metal  ions  are  intro¬ 
duced  into  the  solution  and  can  then  subsequently  hydrolyze  to 
form  hydroxo  complexes.  The  amount  of  metal  ions  obtained  from 
the  solid  can  be  described  by 

|  H20  M+3  +  30H"  (Rll) 

or 

M(OH)3(S)^zi±  M+3  +  30H"  (R12) 

with  the  solubility  product  given  by 

Ksp  -  [M+3]  [0H‘]3  (58) 

It  is  also  possible  that  hydroxo  complexes  will  be  formed 
according  to 

M(OH)3^=±  M(OH)3’n  +  (3  -  n)OH“ 

However,  the  solubility  products  of  such  reactions  are  much 
smaller  than  for  the  former  so  that  these  reactions  will  be 
neglected. 

In  order  to  calculate  the  concentrations  of  the  complexes 
in  solution,  Equation  58  is  solved  in  conjunction  with  Equation  53. 
Therefore,  substituting  Equation  58  for  Cj^  into  55,  the  concen¬ 
trations  of  the  complexes  are 

k 

K_  tt  K, 

k-0,1 . 6  (59) 

The  surface  charge  is  then  determined  as  before  by  the  Equation  15 
for  surface  charging  due  to  multi-component  single  layer  adsorp¬ 
tion. 
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6.3  Surface  Charge  Due  to  Dissociation  of  Surface  Groups 

Oxides  of  tri-valent  metals  such  as  aluminum  and  iron  have 
the  corundum  lattice  in  which  the  oxygen  ions  are  arranged  in 
hexagonal  packing  with  each  metal  ions  coordinated  by  six  oxygen 
ions  and  each  oxygen  ion  by  four  metal  ions.  As  discussed  earlier* 
at  a  freshly  fractured  surface*  the  coordination  of  some  ions 
will  be  incomplete  so  that  hydroxyl  ions  in  solution  will  react 
with  metal  ions  and  protons  with  the  oxygen  ions.  In  this  manner* 
hydroxyl  groups  are  formed  on  the  solid  surface  and  can  ionize 
as  bases  or  acids  to  produce  a  surface  charge.  The  degree  to 
which  the  surface  is  disturbed  determines  how  hydroxy lated  the 
surface  is.  Those  parts  of  the  surface  severely  distorted  are 
hydroxytaced  completely  to  M(OH)3  while  less  distorted  areas  are 
only  partly  hydroxylated  to  MO*OH  (O'Connor  et  al.,  1956). 

The  extent  to  which  the  surface  hydroxyl  groups  ionize  as 
acids  or  bases  depends  on  how  polarized  they  are  by  the  coordin¬ 
ating  metal  ions.  Because  MO-OH  has  two  bonds  and  M(OH)3  three, 
it  is  expected  that  the  hydroxyl  groups  are  more  polarized  in 
MO*OH  than  in  M(OH)3  suggesting  that  the  character  of  MO-OH  is 
acidic  acid  and  M(OH)3  basic  (O’Connor  et  al.,  1956).  On  this 
basis*  MO*OH  can  dissociate  according  to 

Kpl 

MO+(S)  +  2H20  v^=±MO-OH(S)  +  H30+(aq) 

Kp2  ^Rl3^ 

MO-OH(S)  +  H20  ±  M02(S)  +  H30+(aq) 

and  M(0H)3  according  to 
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s 


+3  cl  ,,  , 

M+J(S)  +  2^0  M(0H)+2(S)  +  H30+(aq) 

kc2 

M(0H)+2(S)  +  2H20^=±  M(0H>2 (S)  +  H30+(aq) 

kc3 

M(0H>2  (S)  +  2H20^=i  M(0H)3(S)  +  H30+(aq) 

Kc4 

M(0H)3(S)  MD"(S)  +  H30(aq) 


(R14) 


where  the  subscripts  p  and  c  denote  the  partially  and  conq>letely 
hydrated  oxide,  respectively. 

The  negative  group  MO^  was  chosen  instead  of  M(0H)^  because 
it  is  unlikely  that  a  coordination  number  of  seven  would  be 
realized. 


From  this  mechanism,  it  is  seen  that  at  low  pH  values  the 
surface  charge  will* be  positive  because  the  backward  reactions 
are  favored  while  at  high  pH  it  will  be  negative  because  the 
forward  ones  are  favored. 

Due  to  electrostatic  interaction,  removal  of  a  proton  becomes 
progressively  more  difficult  as  the  charge  on  the  surface  group 
decreases  from  a  high  positive  value  to  a  negative.  This  suggests 
the  following  sequence  for  the  magnitude  of  the  equilibrium  con¬ 
stants 


V>KP2 

Kcl  >  Kc2  >  Kc3  =■  Ko4 
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Also,  at  the  Isoelectric  point  assuming  only  the  presence  of 
uni- valent  species 


K 


Pi 


K 


P2 


[md-oh]  [h+]  „e,„ 

[S’] 

[MO* OH] 


■ZP,C.?P 


so  that 


[h+] 


2 

2pC,p 


[mo+] 

w  [Sfi 


[MOjJ. 


KplKl 

since  for  zero  surface  charge  |^MO+J 
Similarly, 

[H+Jzpc,c  -  Kc3KcA 

Thus,  the  ratio  of  the  isoelectric  points  is 


KPlKp2 


(60) 


(61) 


4* 

H 

zpc,p 

H+: 

zpc,c 

(62) 


As  indicated  above,  hydroxyl  groups  are  more  polarized  in  M0*0H 
than  in  M(OH)^  which  means  that 

V  >  Kc3 


V  ”  Kc4 


Therefore, 


[«+]zpe,p  >  ["+Jzpc 


(63) 


On  the  basis  of  these  facts,  one  concludes  that  metal  oxides 
which  are  extensively  hydrated  to  form  a  surface  layer  of  pre¬ 
dominately  MCOH)^  have  their  isoelectric  points  in  basic  solution 
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and  are  positively  charged  in  neutral  solution;  on  the  other 
hand,  metal  oxides  which  are  partially  hydrated  have  their  iso¬ 
electric  points  in  acid  solution  and  are  negatively  charged  in 
neutral  solution.  For  intermediate  degrees  of  hydration,  the 
surface  consists  of  a  mixture  of  M(OH)3  and  MO-OH  which  accounts 
for  the  wide  range  of  isoelectric  points  of  oxides  and  its  depen¬ 
dence  on  the  treatment  of  the  material  such  as  grinding,  heating, 
ignition,  and  aging  under  water.  These  conclusions  are  in  agree¬ 
ment  with  experimental  observations  for  alumina  or  A1203  (O'Connor 
et  al.,  1956;  Parks,  1965;  Robinson  et  al.,  1964).  " 

The  magnitude  of  the  surface  charge  can  be  expressed  quan¬ 
titatively  by  proceeding  as  in  the  previous  analyses.  By  noting 
the  similarity  of  the  reaction  mechanisms  for  dissociation  o^E 
surface  groups  and  hydroxo  complex  formation,  the  surface  charge 
due  to  dissociation  can  be  easily  deduced.  The  total  charge  due 
to  the  complexes  is 

7 

q  ■  I  zieci 

i-l 

where  c±  is  given  by  Equation  57.  Multiplying  both  numerator 
and  denominator  of  this  expression  by  Lh+J3  one  obtains 

1W*3 
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In  order  to  eliminate  the  extra  reactions  in  hydroxo  complex 
formation  not  present  in  surface  dissociation,  one  sets 

KL  «  K2  -  oo 

K3  "  Kpl 

K4  -  V 

K5  ’  K6  -  0 

c0  ‘  Np 

for  the  surface  charge  contribution  due  to  dissociation  of  MO'OH, 
and 


K, 


-  K 


cl 


K2  *  Kc2 
K3  “  Kc3 


k5  -  0 

K6  -  0 

c0  -  NC 
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r 


for  the  surface  contribution  due  to  dissociation  of  M(OH)g.  Np 
and  Nq  denote  the  total  number  of  sites  which  'are  partially  and 
completely  hydrated  to  provide  the  available  sites  for  the  dis¬ 
sociation  of  MO'OH  and  M(OH)g,  respectively.  Thus,  the  surface 
charge  due  to  dissociation  is  given  by 


N  e 
P 


k_,  r„+i 


pi 


N  e 
c 


TO^cS 


[H+J  I 

3  [h+]  3  ,  2  1h+]  2  .  [h+J  ,  K4  ' 
"TO  ‘  V  , 

2  [atj  .  ,  .  K4  I 

+  tct+  1  +  ivr 

LH  Jj 


TO1^ 


[h+]  .  U+J 

1  W  1  V  1  +  Tr - T7 


TO 


c3 


(65) 


Due  to  the  fact  that  these  dissociation  reactions  are  occur- 
ing  on  a  charged  surface,  the  values  of  the  equilibrium  constants 

i  , 

will  be  a  function  of  the  surface  charge  because  it  becomes  pro¬ 
gressively  more  difficult  to  remove  a  proton  from  the: surface  as 
the  surface  charge  decreases  from  a  positive  value  to  a  negative 
one.  i 

i 

This  effect  on  the  equilibrium  constants  can  be  determined  ! 
if  it  is  recalled  from  classical  thermodynamics  that  the  equili¬ 
brium  constant  and  free  energy  change  AG  per  molecule  of  a 
reaction  are  related  by  ,,  1 


AG  =  -  kT  In  K 


(66) 


By  determining  the  change  in  AG  as  a  result  of  reaction  pn  a 
charged  surface,  the  corresponding  change  in  K  is  obtained. 

Although  the  surface  charge  consists  of  discrete  charged, 
it  is  mathematically  convenient  to  treat  this  charge  as  being 
uniformly  smeared  over  the  surface  (Tanford»  1961).  This  assump¬ 
tion  means  that  the  free  energy  of1 a  surface, molecule  is  being 
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I 


I 


considered  equal  to  the  ayerage  free  energy  of  the  surface  per 
molecule.  Thus,  when  a  proton  of  charge  ze  is  removed  from  a 
surface  of  charge  q,  the  contribution  to  the  change  dG  in  the 
free  energy  G  of  the  surface  due  to  electrical  interaction  is 


dGel 


‘dG 

3q  ze 


Sine 6  the  free  energy  is  given  by  the  work  to  charge  the  surface, 

'  rq 

Gei  “  J  !  Hq/)dq/  / 


Differentiation  under  the  integral  sign  yields 

I  ; 

,3G  , 

.  '  “So1-  *  *0  ' 


so  that 


-  zet, 


Therefore,  if  AG°  is  the  free  energy  change  per  molecule  of  the 
reaction  ocCuring  on  ah  uncharged  surface,  then  AG  for  a  charged 
surface  of  electrical  potential  <i»Q  is  given  by 

AG  -  AG°  -  ze^Q  (67) 

Combining  Equation  66  and  67,  the  equilibrium  constant  can 
be  written  as  ■  i 


K  ■  exp 


-  AG°  +  zetQ\ 

ET  J 


where  the  presence  Of  the  electrical  surface  potential  reflects 
the  surface  chafge  dependence.  The  relationship  between  the 
surface  charge  and  surface  potential  has  been  discussed  in 
Section  4.1. 
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T.f  AG°  is  replaced  by  its  corresponding  equilibrium  constant 
for  reaction  on  an  uncharged  surface,  then 


K  -  Kq  exp  (-j^) 

where 

Ko  -  exp  (-  |f-) 


(69) 
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7.  EXPLANATION  of  the  surface  charge  versus  concentration  curve 


One  Is  finally  in  a  position  to  explain  the  dependence  of 
the  adsorption  density  of  potential  determining  ions  and  the 
resulting  surface  charge  on  these  ions  and  the  ionic  strength. 


When  the  surface  charge  is  a  result  of  adsorption  from 
solution,  it  is  seen  from  the  multi-component  adsorption  isotherm 
given  by  Equation  15  that  the  charge  varies  with 


♦i  +  zic,ol 

- ET - J 


When  the  surface  is  positively  charged,  positive  ions  are  being 
adsorbed  so  that  z^  and  Jq  are  positive  and  thus  also  z^e^Q. 

Since  the  specific  adsorption  energy  is  negative,  a  decrease 
in  makes  AG^  more  negative  and  thus  increases  b^  and  accord¬ 
ingly  also  the  surface  charge.  Since  an  increase  in  ionic  strength 
decreases  the  double  layer  thickness  and  thus  also  the  surface 
potential,  the  above  reasoning  shows  that  the  surface  charge 
increases  with  increasing  ionic  strength  as  found  experimentally. 
Similarly,  for  a  negatively  charge  surface,  the  term  z^etg  is 
again  positive  which  shows  that  the  surface  charge  becomes  more 
negative  with  increasing  ionic  strength.  Also,  since  -  0 
at  the  isoelectric  point,  it  is  seen  that  this  point  is  indepen¬ 
dent  of  ionic  strength  in  agreement  with  experimental  observation. 

In  the  case  when  surface  dissociation  determines  the  surface 
charge,  one  notes  that  the  equilibrium  constant  is  a  function  of 
the  surface  potential  according  to 


K  -  Ko«xp(^0) 

where  z  -  1  for  H+. 

For  a  positively  charged  surface,  an  increase  in  ionic 
strength  makes  |q  less  positive  which  in  turn  decreases  K.  This 
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means  that  the  backward  acid-base  reactions  occur ing  on  the  sur¬ 
face  are  favored  so  that  the  charge  becomes  more  positive. 
Similarly,  for  a  negatively  charged  surface  increasing  ionic 
strength  makes  ♦q  less  negative  which  in  turn  increases  K.  This 
means  that  the  forward  reactions  are  favored  so  that  the  charge 
becomes  more  negative.  Also,  at  the  zero  point  charge,  the  sur¬ 
face  potential  is  zero  so  that  the  equilibrium  constants  and 
thus  the  isoelectric  point  are  independent  of  ionic  strength. 
Again,  these  conclusions  are  in  agreement  with  experimental  obser 
vations . 
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8.  VERIFICATION  OF  THEORETICAL  MODELS 


8 . I  Comparison  of  Adsorption  Theory  with  Experimental  Data 

The  theory  presented  for  the  adsorption  of  ions  from  solution 
onto  solid  surfaces  will  now  be  compared  with  experimental  data. 
The  studies  of  Ottewill  and  Rastogi  (1960)  will  be  used  for  this 
purpose.  They  measured  the  zeta  potential  of  silver  iodide 
solutions  at  pi  4  in  the  presence  of  several  cationic  poly¬ 
electrolytes.  The  radius  of  the  particles  was  200  A. 

In  the  absence  of  polyelectrolytes,  the  particles  had  a 
value  of  ^ a  m  0.66  and  a  zeta  potential  of  -80  mv  due  to  excess 
iodide  ions  in  their  surfaces.  It  is  assumed  that  the  contri¬ 
bution  to  the  surface  charge  due  to  this  excess  of  surface  iodide 
ions  is  not  significantly  affected  by  the  addition  of  a  poly¬ 
electrolyte.  Thus,  the  system  was  treated  as  the  adsorption  of 
one  species  onto  particles  with  an  initial  given  charge. 

The  ionic  strength  was  calculated  from  the  concentrations 
of  Agl  and  the  polyelectrolyte  in  solution.  The  surface  charge 
and  zeta  potential  were  related  through  Equation  35  for  a  planar 
surface.  Although  the  value  of  *a  does  not  warrant  the  use  of 
Equation  35,  it  was  nevertheless  used  because  it  became  too  time- 
consuming  to  utilize  the  computer  program  for  a  spherical  surface. 
This  was  a  result  of  the  large  number  of  calculations  required 
for  use  as  input  to  the  regression  analysis.  The  adsorption 
isotherm  used  for  predicting  the  surface  charge  was  that  given  by 
the  multi-layer  adsorption  Equation  22  for  a  single  component. 

This  isotherm  was  used  because  it  is  known  that  polyelectrolyte 
molecules  can  adsorb  onto  themselves  to  form  a  second  adsorption 
layer.  Also,  the  total  adsorption  energy  term  was  assumed  to  be 
just  the  specific  adsorption  energy  with  the  exclusion  of  the 
surface  potential  term  because  the  latter  is  small  in  comparison 
with  the  former  term.  Under  these  assumptions,  the  model  simpli¬ 
fies  to  one  similar  to  that  proposed  in  the  papers  of  Ottewill 
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which  is  limited  to  C  <  25  mv  and  one  component  system. 

The  values  of  the  parameters  Nq,  and  AGj^  Were 

obtained  by  regression  analysis  subject  to  the  constraint  that 
the  C  versus  c  curves  passed  through  the  measured  isoelectric 
point.  A  computer  program  was  developed  for  this  purpose.  The 
values  obtained  for  the  available  number  of  surface  sites  and  the 
adsorption  energies  were  on  the  order  of  lO1^  sites/cm^  and 

j 

10  Kcal/mole  which  is  in  agreement  with  commonly  accepted  values. 
Also,  the  adsorption  energies  for  the  first  layer  were  always 
greater  in  absolute  magnitude  than  those  for  the  second  which  is 
what  one  would  expect. 

The  C  versus  c  curves  as  predicted  by  the  theory  is  compared 
with  experimental  data  in  Figures  12  to  17  for  various  poly¬ 
electrolytes.  It  is  seen  that  there  is  in  general  good  agreement 
between  theory  and  experiment.  The  discrepancies  which  do  occur 
in  some  of  the  figures  are  around  10“^  M  where  C  is  overestimated. 
This  is  a  result  of  treating  the  surface  as  planar  and  to  the 
lesser  neglecting  the  surface  potential  contribution  to  the 
adsorption  energy  and  can  be  shown  qualitatively  by  using  the 
Debye-Huckel  approximation  (Equation  41).  The  potential  due  to 
a  surface  charge  for  a  spherical  surface  is 

*0S  c(l  +Va) 
and  for  a  planar  one 


The  ratio  of  these  two  equations  gives 

*0S  "  1  +  Ka  *0P  n° ) 

Since  the  Agl  concentration  was  maintained  at  10’^  M  which 
corresponds  to  *  a  ■  0.66,  it  is  seen  that  *a  does  not  vary  signi¬ 
ficantly  until  the  polyelectrolyte  concentration  approaches 
10  M.  In  this  concentration  range,  the  *a  term  acts  like  a 

1 1 T  RESEARCH  INSTITUTE 


61 


63 


ZETA  POTENTIAL  VERSUS  CONCENTRATION  FOR  COLLOIDAL  Agl  IN  THE 
PRESENCE  OF  OCTYLPYRIDINIUM  BROMIDE  (OpyB) 


64 


Figure  14 

ZETA  POTENTIAL  VERSUS  CONCENTRATION  FOR  COLLOIDAL  Agl  IN  THE 
PRESENCE  OF  DODECYLTRIMETHYL- AMMONIUM  BROMIDE  (DTAB) 


Theoretical  Curve: 

=  -8523  cal/raole, 
4G^j®  =  -6149  cal/mole, 
Nn  =  2.005x10^  sites/co 
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-10,717  cal/mole 
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Figure  17 

ZETA  POTENTIAL  VERSUS  CONCENTRATION  FOR  COLLOIDAL  Agl  IN 
PRESENCE  OF  CETYLPYRIDINIUM  BROMIDE  (CPyB) 


constant  and  thus  is  adsorbed  into  the  parameters  during 
regression.  At  polyelectrolytic  concentration  above  10  M, 

^Qg  and  4'qp  become  identical  because  *a  >:>  1*  Therefore,  good 
agreement  between  theory  and  data  is  expected  in  the  concentration 

range  10”^  M  >  c  >  10”^  M  as  found.  However,  in  the  range 

10~3  M  <  c  <  lO”0  M,  Ka  will  vary  and  also  be  on  the  order  of 

unity.  As  a  result,  it  is  seen  from  Equation  70  that  the  assump¬ 

tion  of  a  planar  surface  will  overestimate  the  true  surface 
potential. 

Therefore,  it  may  be  concluded  that  the  theory  presented 
for  the  development  of  surface  charge  due  to  adsorption  is  both 
conceptually  and  quantitatively  correct. 

8.2  Comparison  of  Surface  Dissociation  Theory  with 
Experimental  Data 

The  theory  presented  for  the  development  of  surface  charge 
due  to  the  dissociation  of  surface  groups  will  now  be  compared 
with  experimental  data.  The  study  of  Modi  and  Fuerstenau  (1957) 
will  be  used  for  this  purpose.  They  measured  the  zeta  potential 
of  alumina  in  aqueous  solutions  at  various  pH  values  and  salt 
concentrations  by  the  streaming  potential  method  in  a  packed  bed. 
The  size  of  the  particles  was  48/65  mesh. 

Since  the  particle  size  is  on  the  order  of  microns,  the 
value  of  Ka  is  much  greater  than  one  so  that  the  particles  may 
be  treated  as  planar  surfaces.  Thus,  the  surface  charge  was 
calculated  from  the  zeta  potential  measurements  according  to 
Equation  35.  Also,  since  the  zero  point  charge  occurs  at  pH 
9.4,  it  was  assumed  that  the  surface  consisted  only  of  Al(OH)^ 
groups . 

As  a  result,  the  adsorption  isotherm  given  by  Equation  65 
with  Np  -  0  was  used  to  predict  the  surface  charge.  The  values 
of  the  parameters  Nq,  Kq^,  Kq£»  Kq3»  were  obtained  by 
regression  analysis  performed  by  the  computer  program  to  give 
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I 


No  - 

1.557xl013  sites/cm2 

K01  - 

i.o  •; 

o* 

N5 

1 

0.02138 

K03  “ 

4. 678x10" 6 

o 

-P- 

I 

} 

3. 613x10“ 11 

It  is  noted  that  the  number  of  sites  is  in  agreement  with  the 

!  !  1  O  o 

common  fact  that  it  should  be  on  the  order  of  10  sites/cm  .  1 

l  f 

Also,  the  values  of  the  equilibrium  constants  decrease  in  the 
order  Kq^, . ,Kq^  as  predicted  and  are  of  reasonable  magni¬ 

tudes  .  It  was  found  that  calculations  were  insensitive  to 
which  combined  with  its  large  value  indicates  that  Al+3  was  not 
formed  on  the  surface.  This  makes  sense  because  it  is  not 
expected  that  the  bare  metal  would  be  exposed.  , 

The  surface  charge  predicted  by  Equation  65  as  a  function 

|  » 

of  pH  and  ionic  strength  is  compared  with  experimental  data  in 
Figure  18.  Similar  comparisons  are  given  in  Figure  19  where 
surface  charge  is  plotted  against  pH  and  in  Figure  20  where  sur¬ 
face  charge  is  plotted  against  salt  concentration  for  various 
pH  values.  From  these  graphs,  it  is  seen  that  the  theoretical 
and  experimental  results  are  in  good  agreement  at  high  and  low  t 
pH  values  but  that  there  is  somewhat  pobr  agreement  aroung  pH  7. 
Also,  the  isoelectric  point  is  incorrectly  predicted  to  be  i 
pH  7.8  instead  of  the  experimental  value  of  pH  9.4.  However, 
these  discrepancies  can  probably  be  accounted  for  by  hypothesizing 
that  a  small  fraction  of  the  AlCOH)^  surface  groups  are  really 
A10*0H  groups.  It  is  expected  that  such  a  surface  would  allow 
the  calculated  isoelectric  point  to  occur  at  a  higher  pH.  This 
projection  is  based  on  numerical  computations  performed  during 
the  course  of  the  regression  analysis. 

i 
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Surface-  Charge.  .  eoul/cm‘ 


n  o 

1  1  Theoretical  curve:  NQ  *  1.557x10  sites/cm  , 
Kq^  -  1.0,  Kq2  -  0.02138,  KQ3  -  4.678xl0‘6, 
Kq4  -  3.613X10"11 

O  0.1  M  •  =0.01  M 

I 

Data : 

P  0.001  M  A  0.0001  M 


CN 


Figure  18 

SURFACE  CHARGE  OF  ALUMINA  VERSUS  pH  AND  IONIC  STRENGTH 

1  i  i 

i 
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Surface  Charge,  ..  coul/cm 


Figure  19 

SURFACE  CHARGE  OF  ALUMINA  VERSUS  pH 
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_  Theoretical  curve: 

N0  -  1.557xl013  sites/cm2, 

*0!  -  1.0,  Kqj  -  0.02138, 

*03  -  4.678x10“®,  -  3.613xlO-11 

O  Experimental  Data 

o 


Figure  20 


Although  there  is  not  complete  agreement  between  theory  and 
data,  there  is  sufficient  agreement  to  substantiate  the  concept 
of  surface  dissociation  as  a  charging  mechanism  for  alumina  and 
the  basic  mathematical  model. 
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9.  TOTAL  ENERGY  OF  INTERACTION  CURVES 


9.1  van  der  Waals-London  Attraction 

As  discussed  in  che  literature  review,  the  total  interaction 
energy  of  two  colloidal  particles  is  given  by  the  sum  of  the 
Interaction  energies  due  to  van  der  Waals-London  attractive 
forces  and  to  electrical  double- layer  Interaction. 

The  energy  of  interaction  between  two  spherical  particles 
was  considered  in  detail  in  this  review,  and  it  was  shown  that 
for  small  separation  distances 

*  Afl i  Ao 

va  *  +  v2m0  (71) 

where 

A  is  the  Hamaker  constant 
a^,a2  are  the  radii  of  the  particles 

Hq  is  the  shortest  distance  between  their 
surfaces . 

The  effects  of  retardation  and  adsorbed  layers  were  also  con¬ 
sidered. 

The  theory  for  the  van  der  Waals  force  is  well  established 
and  complete  so  that  no  improvements  or  modifications  were 
necessary.  For  this  reason,  there  will  be  no  further  discussion 
of  this  subject,  and  the  reader  is  referred  to  the  literature 
review  for  details.  However,  in  the  case  of  electrical  double¬ 
layer  interaction,  it  was  necessary  to  extend  the  present  theory 
as  given  on  the  following  pages. 

9.2  Electrical  Double-Layer  Repulsion 

Although  there  have  been  detailed  calculations  carried  out 
by  several  investigators  (Kruyt,  1952;  Devereux  and  de  Bruyn, 

1963;  Honig  and  Mul,  1971)  of  the  energy  of  interaction  between 
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Che  electrical  double  layers  associated  with  two  particles,  the 
results  of  these  studies  are  given  in  tabular  form  rather  than  as 
equations  because  graphical  or  numerical  integrations  were  per* 
formed.  However,  since  this  interaction  energy  is  to  be  used 
later  in  this  project  in  calculations  which  involve  integration 
and  the  solution  of  differential  equations,  such  solutions  will 
make  these  further  operations  extremely  unwielding.  Thus,  there 
is  a  distinct  need  for  equations  which  can  approximately  describe 
this  Interaction.  Unfortunately,  the  two  moat  promising  ones 
available  are  restricted  to  the  interaction  of  particles  with 
the  same  surface  potential  (Honig  and  Hul,  1971)  whi^e  for 
particles  with  different  potentials  to  the  case  where  the  poten¬ 
tials  do  not  exceed  an  absolute  value  of  50  mv  (Hogg  et  al., 

1966).  In  this  section,  the  limits  of  these  approximate  equations 
will  be  extended. 

9.2.1  Approximation  for  Large  and  Moderate  Particle 
Separation 

The  basis  of  this  method  was  originally  worked  out  by  Verwey 
and  Overbeek  (1948)  for  two  fiat  plates  with  identical  surface 
potentials.  When  the  two  surfaces  are  relatively  far  apart, 
there  is  small  interaction  so  that  the  potential  between  them 
can  be  obtained  by  adding  the  potentials  due  to  two  isolated 
surfaces.  From  this,  the  interaction  between  spheres  is  then 
obtained  by  Derjaguin's  method  which  considers  the  spherical 
surface  to  be  made  up  of  flat  rings.  Justification  for  this 
procedure  *s  that  for  the  case  of  identical  surface  potentials, 
results  of  this  method  are  in  good  agreement  with  exact  numerical 
calculations  (Honig  and  Mul,  1971). 

An  expression  for  the  potential  of  an  isolated  flat  plate 
in  obtained  from  the  Poisson  Boltzmann  equation  (Equation  31) 
written  in  the  form  (Verwey  and  Overbeek,  1948) 
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(72) 


•  Sinhy 
d5z 


where 


y 

5 


subject  to 


zet 

XT 

<X 


y0  st  5  -  0 


and 


-  d 


S3f 


0  as  5  -® 


The  hyperbolic  sine  term  is  a  result  of  assuming  that  nil  the 
electrolytes  in  solution  are  of  valency  z.  This  slmpllcatlon  is 
permissible  because  ions  with  the  same  charge  of  the  surface  are 
not  important  (Kruyt,  1952). 

The  solution  of  Equation  72  is 


ey/2  -  +  l  ±  -  »«'S 

.y«“  ♦  i  -  (.y°,Z  -  1 


(73) 


Writing  Equation  73  as 


2  In 


1  +  ve”3 


1  -  ve 


(74) 


where 


v  - 


cy°/2  -  1 

7^77 


it  is  seen  that  sufficiently  far  away  from  the  surface  such  that 
5  »  1  Equation  74  simplifies  to 

y  ~  2  ln(l  +  2ve’^) 


4ve 


-5 


(75) 
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as  a  result  of  expansion  in  series  of  the  logarithmic  term.  It 
has  been  shown  that  Equation  75  la  a  very  good  approximation  for 
5  >  1  (Verwey  and  Overbeek,  1948). 

If  the  two  plates  are  located  a  distance  2d  apart »  then 
from  Equation  75  the  contribution  of  Plate  1  at  x  -  0  to  the 
potential  in  the  solution  is 

yl  “  4vle’^  (76) 

and  that  of  Plate  2  at  x  -  2d 


y2 


4v2*«-2‘<» 


on 


Thus,  adding  Equations  76  and  77,  the  potential  between  the 
plates  is 

y  -  +  Vje  (^“2*  d)  (73) 

The  force  and  thereby  the  interaction  energy  between  the 
plates  are  obtained  by  following  the  method  by  Verwey  and  Over¬ 
beek  (1948) .  It  is  assumed  that  the  plates  are  kept  apart  by  a 
pressure  P  and  the  liquid  oetween  them  is  in  contact  with  a  large 
reservoir  at  potential  zero  and  pressure  Pq.  A  force  balance  on 
any  point  in  the  liquid  between  the  plates  gives  the  force  F  at 
that  point  as 


dF  ■  dp  -  pEdx  ■  dp  +  pdt  ■  0 

Introducing  the  Poisson  equation,  Equation  79  becomes 


(79) 


Integration  yields 


constant 


(80) 
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The  constant  of  integration  is  given  by  the  pressure  where  •  0. 
This  point  xm  is  found  by  taking  the  derivative  of  Equation  78 
with  respect  to  ?  and  setting  it  equal  to  zero. 


or 


where 


&  - 


"5, 


m 


+  v^e 


V2<dl 


-  0 


2K  (xm-d) 


m 


m 


-  Kd 


♦* 


In 


v2 


<d  +  6 


0 


(81) 


(82) 


The  force  driving  the  plates  apart  due  to  the  double  inter¬ 
action  is  the  force  relative  to  the  reservoir  pressure  and  is 
obtained  from  Equation  79  as 


♦-0  0 

Writing  the  charge  density  in  accordance  with  Equation  72  as 
P  “  2zenQSinh 


and  substituting  it  into  Equation  83, 

zet 


P  “  2nQkT  cosh  ^  -  1 


(84) 
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-'lnce  weak  interaction  between  the  particles  has  been  assumed, 
it  is  sufficient  to  use  the  first-order  term  in  the  Taylor  expan¬ 
sion  of  Equation  84  about  ♦  »  0  as  given  by 


VUt r 


nokTV 


(85) 


Substitution  of  Equation  78  into  Equation  85  yields 
„  ,2  F  2  “2?m  2  2(5.  -2*d) 

P  -  "okTym2  ■  16  [vj2.  ”'  +  v22e  »  +  2v1v2e-2Kd  J  n0kT 

,kTe’ 


16 


[vj2e- 


“  +  V*2*  ^  2v1v2 


1  ,  ~  -2*d 

J  n0l 


(86) 


The  energy  of  interaction  VR  is  determined  from  the  relation 


2d 

VR  -  -  /  P(2d)d(2d) 

oo 


which  upon  substitution  of  Equation  86  gives 

VR  •  *  16  [vl2**2"  +  V*26  +  2vlv2  ]n0kT  / 

oc 


16 


v  V26 
1  e 


+  V 


^e26  +  2vtv 


ln0kT 

.  - 


-2<d 


(87) 


From  this  result,  it  is  now  possible  to  calculate  the  inter¬ 
action  between  spherical  particles  by  means  of  Derjaguin's  method 
(Verwey  and  Overbeek,  1948;  Kruyt,  1952).  If  the  thickness  of 
the  double  payers  is  small  relative  to  the  radius  of  the  particles 
(*a  >  5),  the  double  layer  interaction  between  spherical  particles 
can  be  considered  to  consist  of  contributions  from  differential 
parallel  rings  as  shown  in  Figure  21.  Treating  each  pair  of  rings 
independently  as  two  flat  plates, the  energy  of  interaction  between 
two  spherical  particles  can  be  written  as 
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Figure  21 

GEOMETRICAL  CONSTRUCTION  USED  IN  THE  CALCULATION 
OF  THE  INTERACTION  BETWEEN  TWO  DISSIMILAR  SPHERICAL 
PARTICLES,  RADII  AND  a2  FROM  THE  INTERACTION 

OF  TOO  INFINITE  FLAT  PLATES 
(from  Hogg  et  al.,  1966) 


a c 


VR  -  / 


0 


2nhvputedh 


(88) 


From  the  geometry  of  Figure  21 

r— 

H  -  a,  +  -  fa^  -  h 


1  ’  M2  "  \/al  '  "2  *  ~\l8 2  *  h2  +  Hq 


Differentiation  of  Equation  89  with  respect  to  h  yields 


dH 


h2/al2  *2  fl*  I  >2/»:2 


hdh 


which  simplifies  to 


hdh  ~ 


ala2 
al  ^  a2 


dH 


(89) 


(90) 
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under  the  reasonable  assumption  h  «  a^  and  h  «  &2 
Substituting  Equation  90  into  Equation  88 


vr  '  r-Hr  1  v 

1  ^  Hrt 


Plate 


(H)dH 


(91) 


Finally,  substituting  Equation  87  into  Equation  91  and  carrying 
out  the  integration  with  H  -  2d 

nnkT  -*Ha 


32rra]®2  T  7  - 

VR  ’  T^TT2  L  V' 


2*  +  *2V4  +  2 


]  V 

vi  v2  J  "TJ 


(92) 


For  the  special  case  where  ai  “  a2  "  a>  Equation  92  reduces  to 


64TtavZn0kT  -KHQ 

- 77 -  e 


(93) 


which  has  been  extensively  compared  by  Honig  and  Mul  (1971)  and 
found  to  be  in  good  agreement  with  exact  numerical  calculations 
for  kHq  >  1.  Since  van  der  Waal's  forces  start  to  dominate  around 
*Hq  •  1,  Equation  92  is  a  very  useful  approximation. 

^•2.2  Approximation  for  Small  Particle  Separation 

The  rationale  for  this  method  is  that  for  small  separation 
distances  between  two  particles  the  potential  in  solution  will 
not  decrease  a  great  amount  from  that  at  the  surface.  Thus,  one 
can  linearize  the  Poisson  Boltzmann  equation  by  expanding  the 
charge  density  about  the  value  of  the  surface  potential.  The 
procedure  from  this  point  on  is  similar  to  that  used  in  obtaining 
the  previous  approximation.  The  interaction  between  two  flat 
plates  is  found  first,  and  then  Derjaguin's  method  is  used  to 
calculate  the  interaction  between  two  spherical  particles.  Hogg, 
Healy,  and  Fuerstenau  (1966)  have  carried  out  this  method  for  the 
case  where  the  charge  density  was  expanded  about  a  potential  of 
zero.  Their  results  simplified  the  calculations  for  the  present 
task. 


1 1 T  RESEARCH  INSTITUTE 


81 


It  has  been  commonly  assumed  in  the  literature  for  mathe¬ 
matical  convenience  that  as  two  particles  approach  each  other 
their  surface  potentials  do  not  change*  and  as  a  result*  their 
surface  charges  diminish.  This  implies  that  adsorption  equili¬ 
brium  is  rabidly  established.  However,  it  is  doubtful  that  this 
is  true,  and  a  more  appropriate  assumption  is  constancy  of  the 
surface  charge  although  the  actual  case  lies  somewheie  in  between 
these  limiting  cases.  This  implies  slow  desorption  which  seems 
more  reasonable  in  light  of  the  short  particle  contact  times 
involved  due  to  Brownian  motion.  Unfortunately,  this  assumption 
complicates  the  mathematics  somewhat.  For  this  reason,  the  fact 
that  the  results  of  either  assumption  do  not  differ  significantly 
will  be  used  so  that  constant  surface  potentials  can  be  assumed. 
This  has  been  shown  to  be  true  by  Honig  and  Mul  0  971)  who  per¬ 
formed  exact  numerical  calculations  for  the  interaction  energy 
at  both  constant  surface  potential  and  charge.  Only  at  very 
small  separation  distances  where  the  van  der  Waals  forces  would 
dominate  are  the  results  in  disagreement. 

The  first  step  in  the  derivation  is  to  linearize  the  Poisson- 
Boltzmann  equation 


-J^et/kT 

zieniOe 


(94) 


Expanding  in  series  the  exponential  terms  about  some  reference 
potential  such  as  the  average  of  the  potentials  of  the  two 
particles  and  retaining  only  the  linear  terms,  Equation  94  becomes 

V2*  -  -  Pr  +**(*  -  *r)  (95) 

where 
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eniOe 


-ZiCtR/kT 


-z,e*D/kT 


4neA  V  2  a  *ic’R 
cTTT  L  zi  niOe 


Equation  95  is  to  be  solved  for  two  flat  plates  subject  to  the 
boundary  conditions 


♦  -  ♦ 


*  -  ♦ 


01 


02 


at  z  ■  0 

at  x  ■  2d 


If  a  change  of  variables  is  made 
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(96) 


then  the  problem  to  be  solved  is 
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dx 
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(97) 


x  *  0 


x  *  2d 


The  solution  is 


A  A  A 

♦  *  cosh  *x  + 


♦q2  “  ^oi  2^d 

Sinh  2^d 


Sinh*x 


(98) 
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The  surface  charge  is  calculated  from  the  relation 

i 

which  was  derived  earlier  and  yields  the  result 

r A  A 

"  *  7m  [^02  cosec^(2<d)  * 


01  coeh(2Kd)] 

(99) 

qj  cosech(2<d)J 

(100) 

Since  the  surfape  charge  and  potential  are  linearly  related, 
the  free  energy  G  can  be  determined  from 

G  “  -  /  q(t)dt 

1  0 

G  -  -,7[qlf01  +V02] 

with  the  result  that  the  interaction  energy  Vplflte  is 

1 

VPlate  "  G2d  ‘  ^ 

•  W  |[*012  +  *022  -  4*r(*01  +  *02)]  [l  -  coth(2'Jd)] 

+  [2*01*02  •  Mr(*01  +  *02)]  cosech(2Kd)|  (101) 


where 


Mr  -  tR  +  ^7  PR 


1  Then,  using  Derjaguin's  method,  the  interaction  energy  for 
two  spherical  particles  is 


.  ! 
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I 


(102) 


For  fR  ■  0  and  thus  AfR  *0  and  <  ■  <  ,  Equation  102  simplifies 
to  that  obtained  by  Hogg  et  al.  (1966)  which  they  showed  to  be  In 
agreement  with  exact  numerical  results  for  Iq  <  60  mv.  Thus,  It 
Is  expected  that  Equation  102  will  be  a  good  approximation  for  the 
case  where  the  separation  distance  is  such  that  the  potential  in 
solution  does  not  vary  more  than  60  mv  from  the  surface  potential. 
This  situation  should  occur  for  at  least  *Hq  <  1  If  not  a  greater 
value.  The  reasoning  behind  this  is  that  since  the  double  layer 
thickness  of  an  isolated  particle  is  l/< ,  *Hq  -  1  suggests  that 
potentials  at  positions  no  greater  than  half  a  double  layer  thick¬ 
ness  are  being  considered.  At  a  distance  of  half  a  double  layer 
away  from  the  surface,  the  Debye-Huckel  approximation  f  ■  X 

shows  that  the  potential  has  only  diminished  to  60  percent  of 
that  at  the  surface.  Thus,  for  a  maximum  surface  potential  of 
150  mv,  the  potential  only  reduces  to  90  mv  so  that  the  change 
is  within  60  mv.  Therefore,  since  Equation  92  is  valid  for 
<H0  >  1  and  Equation  102  for  *Hq  <  1,  these  two  equations  together 
provide  a  good  representation  for  the  entire  energy  of  inter¬ 
action  curve  for  two  spherical  particles. 
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10.  CONCLUSIONS  AND  SUMMARY 


It  has  been  shown  that  the  Gouy-Stirn-Grahame  model  is 
Inadequate  for  explaining  the  observed  dependence  of  the  zeta 
potential  and  surface  charge  on  ionic  strength,  surfactant  con¬ 
centration  and  pH.  It  was  necessary  to  explain  two  common  types 
of  experimental  data:  (1)  zeta  potential  versus  concentration  and 
(2)  surface  charge  versus  ionic  strength  versus  pH.  It  was 
shown  that  these  phenomena  are  governed  by  competing  phenomena 
(1)  adsorption  of  ions  from  solution,  (2)  dissociation  of  surface 
groups  due  to  chemical  surface  reactions,  and  (3)  the  thickness 
of  the  electrical  double  layer  due  to  counter  ions. 

In  the  case  of  adsorption,  the  following  was  found.  At  low 
concentration,  the  adsorption  effect  dominates,  and  the  zeta 
potential  increases  in  absolute  value  with  increasing  concentra¬ 
tion.  At  high  concentrations,  the  double  layer  thickness  effect 
dominates,  and  the  zeta  potential  approaches  zero  with  increasing 
concentration.  The  zeta  potential  versus  concentration  curve  can 
be  predicted  mathematlcaly  by  using  a  multi-component  adsorption 
isotherm  in  conjunction  with  a  relationship  between  the  zeta 
potential  and  surface  charge.  For  single  layer  adsorption. 
Equation  15  is  used  for  the  isotherm  while  for  multi-layer  adsorp¬ 
tion,  Equation  22  is  used.  For  planar  surfaces,  the  zeta  poten¬ 
tial  and  surface  charge  are  related  through  Equation  35  while  for 
spherical  surfaces,  a  computer  program  serves  this  purpose.  When 
the  species  to  be  adsorbed  undergoes  hydroxo  complex  formation, 
it  is  necessary  to  calculate  the  concentration  of  the  various 
complexes.  For  an  externally  added  electrolyte,  Equation  57  is 
used  while  for  ions  derived  from  the  particle,  Equation  59  is 
used.  Then  Equation  15  is  used  to  compute  the  surface  charge  due 
to  adsorption  of  the  complexes.  Predicted  zeta  potential  versus 
concentration  curves  were  shown  to  be  in  good  agreement  with 
experimental  data  in  the  literature  for  Agl  sols  in  the  presence 
of  pc l /electrolytes . 
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For  the  case  of  surface  dissociation,  the  following  was 
found.  For  metal  oxides,  the  surface  consists  of  hydroxides  of 
varying  degrees  of  hydration.  The  greater  degree  of  hydration 
the  higher  ia  the  pH  at  which  the  zero  point  charge  occurs. 
Equation  65  is  used  to  calculate  the  surface  charge  while 
Equation  69  gives  the  dependence  of  the  equilibrium  constants  on 
the  surface  charge  or  zeta  potential.  The  zeta  potential  and 
surface  are  related  as  described  above.  It  was  shown  that  this 
dependence  of  the  equilibrium  constants  is  the  reason  why  the 
surface  charge  increases  in  absolute  magnitude  with  increasing 
ionic  strength.  Predicted  surface  charge  versus  ionic  strength 
versus  pH  curves  were  shown  to  be  in  good  agreement  with  experi¬ 
mental  data  in  the  literature  for  alumina  at  various  pH  and  NaCl 
solutions . 

In  order  to  predict  the  interaction  of  two  particles  when 
they  approach  each  other,  the  theory  of  interaction  energy  curves 
was  studied.  It  was  found  that  the  theory  for  van  der  Waals 
forces  is  sufficiently  complete  so  that  no  improvements  were 
necessary.  However,  in  the  case  of  the  electrical  forces  due  to 
double  layer  interaction,  it  was  necessary  to  extend  the  present 
theory.  Although  exact  numerical  calculations  have  been  per¬ 
formed,  they  are  not  c  '•nvenient  to  use  so  that  it  was  necessary 
to  derive  approximate  analytical  expressions.  The  entire  inter¬ 
action  energy  curve  for  spherical  particles  of  arbitrary  size 
and  surface  potential  is  described  by  Equation  102  for  small 
separation  distances  (*a  <  1)  and  Equation  92  for  moderate  and 
large  separation  distances.  These  equations  are  in  good  agree¬ 
ment  with  the  limited  exact  results . 

When  the  attraction  energy  curve  due  to  van  der  Waals-London 
forces  and  the  repulsion  energy  curve  due  to  double  layer  inter¬ 
action  are  added  together,  the  total  interaction  energy  curve 
is  obtained  from  which  the  stability  of  a  colloidal  system  can 
be  deduced. 
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